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A  uniformly  differentiable  approximation  scheme 
for  delay  systems  using  splines 

K.  iTO^and  F.  Kappel^ 

Center  for  Control  Sciences,  Brown  University 

and 

Institute  for  Mathematics,  University  of  Graz 

Abstract.  A  new  spline-based  scheme  is  develop>ed  for  linear  retarded  functional  differential 
equations  within  the  framework  of  semigroups  on  the  Hilbert  space  R”  x  .  The  approxi¬ 
mating  semigroups  inherit  in  a  uniform  way  the  characterization  for  differentiable  semigroups 
from  the  solution  semigroup  of  the  delay  system  (e.g.  among  other  things  the  logarithmic 
sectorial  property  for  the  spectrum).  We  prove  convergence  of  the  scheme  in  the  state  spaces 
R”  X  L-  and  /f’ .  The  uniform  differentiability  of  the  approximating  semigroups  enables  us 
to  establish  error  estimates  including  quadratic  convergence  for  certain  classes  of  initial  data. 
We  also  apply  the  scheme  for  computing  the  feedback  solutions  to  linear  quadratic  optimal 
control  problems. 


1.  Introduction  and  preliminaries 
Ill  this  paper  we  consider  the  linear  hereditary  control  system 


(1.1) 


( 

=  Y^AiT{i  +  B,)  + 

1  =  0 


A{B)T(t  +  6)d9 


x(0)  =  T].  x{6)  =  ip{B)  a.e.  on  (— r,0), 

y(t)  =  C.r(/).  i  >  0. 


t  >  0. 


where  — r  =  <  •  •  ■  <  0o  =  0,  x{t)  €  R",  u(f)  G  R"'  and  y{t)  G  R'’.  Furthermore.  .■!(•) 

is  an  n  X  n  matrix  valued  square  integrable  function  on  [— r, 0].  It  is  well-known  that  for 
(y-r)  £  Z  =  R"  X  £■(  —  /'. 0;  R" )  and  u  G  ex:  R"' )  system  (1.1)  admits  a  unique 

solution  X  G  I^(— r.TiR")  fl  //'(O.TiR")  for  any  T  >  0.  We  define  the  operators  S(i). 
t  >  0,  by 

S(f)(77,vt)  =  (x(t),x,),  t  >  0,  G  Z, 

where  x{t)  is  the  solution  of  (1.1)  with  initial  data  and  u{t)  =  0.  The  segment  j, 

is  defined  by  Xi{6)  =  x{t  -h  B),  —r  <  B  <  0.  The  family  S(-)  is  a  strongly  continuous 
semigroup  on  Z  with  infinitesimal  generator  A  given  by  (see  for  instance  [6]  or  [25]) 


dom  A  =  {(77,(,?)  e  Z  I  ip  e  H^(-r,0;R"),  rj  =  1,5(0)}, 

A(<p(0),i,p)  =  {Ly>,ip)  for  (v?(0),<,5)  G  dom^I, 

'  tVork  done  by  this  author  was  supported  by  AFOSR  under  Contract  No.  F-49G20-8G-C-01  U^y  N.\S.\ 
under  Grant  No.  NAG-1-517  and  by  NSF  under  Grant  No.  UINT-8521208. 

'Work  done  by  this  author  was  supported  by  AFOSR  under  Gram  No.  84-0398  and  by  FWF(Auslria) 
under  Grants  S320G  and  PG005. 
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where  for  continuous  v? 


1=0 

If  we  define  the  input  operator  B :  R""  — +  Z  and  the  output  operator  C:  Z  —*  by 

Bu  =  (Rt/,0),  u  G  R”', 

C{Tj,^)  =  Cr),  (77,9)  ez, 

then  (1.1)  is  equivalent  to  the  following  abstract  system  in  Z: 


Aie)^ie)de. 


z(f}  =  Az(t)  +  Bu(t),  t>0,  ziO)  = 
y(t)=Cz(i),  i>0. 

More  precisely,  a  function  t:  [0,oo)  — »  R"  is  a  solution  of  (1.1)  if  and  only  if  the  function 
z(t)  =  {x{t),xt),  <  >  0,  is  a  mild  solution  of  (1.3),  i.e., 


(1.4)  r(t)  =  5(t)(7/,v?)  +  /  S{i  -  s)Bii{s)ds,  f  >  0. 

Jo 

We  shall  frequently  make  use  of  the  following  facts:  The  spectrum  of  A  is  only  point 
spectrum  and  A  G  cr(A]  if  and  only  if  det  A(A)  =  0,  where  (/„  is  the  ??  x  ??  identity*  matrix') 

A(A)  =  AJ„ -I(e^  /„),  AgC. 

The  resolvent  of  A  is  given  by  {p{A)  denotes  the  resolvent  set  of  A) 

(A/->1)-'(/7,^)  =  (vH0),v'').  AGp(,4). 

where 


(1.5) 


fO 

iie)  =  e^UiO)+  /  e^<^-">9(s)ds,  -r<^<0, 

Je 

m  =  A(A)-'  (77  +  L{J  t^^-^^^{s)ds)). 


Together  with  the  solution  semigroup  S(  )  w'e  shall  need  the  adjoint  semigroup  S(  )*. 
Its  infinitesimal  generator  A*  is  given  by  (see  for  instance  [6,25]) 

^  dom-4*  =  {(y,t'’)  €  Z  ]  w  G  i/^(-r,0;R")  and  =  aJ y] , 

A*{y.x!')  =  (t/iO)  +  Ajy.A'^(-)y-  lij  for  (y,t/')  G  dom>l*, 

where  w  =  V’  +  Ajyx^^r.e,)^  \M  denoting  the  characteristic  function  of  a  set 
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We  shall  also  use  the  state  space  r,0; R")  for  systenr  (1.1).  Endowed  with 

the  inner  product 

=  (s^(0),t/’(0))Rn  + 

is  isomorphic  to  dom^  with  the  graph  norm.  The  isomorphism  i:  dom>l  \s 

given  by 


/(9(0),(^)  =  (^  for  (9(0),9)  G  domA 

for 

Since  S(t),  <  >  0,  restricted  to  domw4  forms  a  Co-semigroup  on  dom^  (with  the  graph 
norm), 

S{t)  =  iS{t)L~^ ,  t  >  0, 

defines  a  Co-semigroup  on  In  fact,  =  i,,  where  x(/)  is  the  solution  of  (1.1)  with 

u  =  0  and  initial  data  (<,5(0),(,?)  G  dom>l. 

If  we  observe  that  /JS(t  —  s)Bu{s)ds  =  (x(t),ij)  for  t  >  0  and  u  G  £j^Qp(0,  oo;  R™ ), 
where  x(/)  is  the  solution  of  (1.1)  with  initial  data  r;  =  0,  (,?  =  0,  then  it  is  not  difficult  to 
see  that  ^ 

t  — >  /  S{i  —  $)Bu{s)  ds 
Jo 

defines  a  continuous  map  into  dom^  with 


/  S(t  -  s)ei/(s)ds|^,  <  ^{t)W\mo,f,K- 
Jo 


for  t  >  0,  where  m;  R"*"  — >  R"*"  is  a  nondecreasing  function.  Hence  (1.1)  with  ?;  =  >r(0). 
G  is  also  well-posed  in  .  In  fact. 


(1-S) 


X,  =  i:{t)  =  S{t)^  -f  / 


fs(t- 

Jo 


s)Buis)  ds. 


t  >  0. 


where  z{i)  is  given  by  (1.4)  and  x(t)  is  the  solution  of  (1.1)  with  initial  data  (;,?(0),!,?)  € 
dom^. 

Approximation  of  hereditary  control  systems  by  ODE-systems  has  some  history  already. 
In  [2]  the  approximating  systems  were  obtained  by  using  the  so-called  averaging  projec¬ 
tions,  i.e.,  projections  onto  a  subspace  of  step-functions.  The  proof  of  convergence  as  in 
almost  all  of  the  following  papers  was  based  on  a  version  of  the  Trotter-Kato  theorem  in 
semigroup  theory.  In  order  to  get  higher  rates  of  convergence  in  [3]  a  scheme  was  devel¬ 
oped  which  uses  first  order  or  cubic  splines  as  approximating  elements.  If  one  considers 
minimization  of  a  quadratic  cost  criterion  for  system  (1.1)  a  very  desirable  property  of  an 
approximation  scheme  is  that  also  the  adjoints  of  the  approximating  semigroups  converge 
strongly  to  the  adjoint  of  the  solution  semigroup  of  the  delay  system  [7].  The  A\’-scheme 
developed  in  [2]  has  this  property  whereais  the  spline  scheme  developed  in  [3]  docs  not  have 
it.  Therefore  in  [15]  a  different  spline  scheme  was  presented  where  also  strong  convergence 
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of  the  adjoint  semigroups  is  guaranteed.  In  contrast  to  the  AV-schcme  [22]  the  spline 
scheme  of  [15]  does  not  have  the  property  that  exponential  stability  of  the  delay  systems 
guarantees  uniform  (with  respect  to  the  approximation  parameter)  exponential  stability 
of  the  approximating  systems.  This  spline  scheme  still  has  the  property  of  uniform  output 
stability  [16]  which  is  enough  in  order  to  deal  with  the  infinite  time  horizon  problem  for 
(1.1)  [17].  In  [12,13,10]  a  scheme  using  Legendre  polynomials  as  approximating  elements 
was  developed.  The  construction  is  based  on  Lemezos’  tau-method.  The  scheme  has  all 
the  qualitative  properties  which  are  stated  above  to  be  true  for  the  AV-scheme.  Also  using 
Legendre  polynomials  a  different  scheme  was  developed  in  [14]  using  the  basic  ideas  of 
[15].  Numerical  evidence  indicates  that  this  Legendre-scheme  has  analogous  properties  as 
the  Legendre-tau-scheme  of  [12],  but  uniform  exponential  stability  for  the  approximating 
systems  has  not  been  established  yet.  Using  piecewise  linear  functions  in  [21]  a  scheme 
was  developed  in  the  spirit  of  [15]  which  also  has  all  the  qualitative  properties  mentioned 
above  for  the  AV  and  Legendre- tau  schemes.  In  all  the  papers  mentioned  up  to  now  no 
convergence  rates  have  been  established  or  only  convergence  rates  which  obviously  are 
not  optimal.  For  instance  in  [3]  only  linear  convergence  for  smooth  data  in  case  of  first 
order  splines  is  established.  In  a  recent  paper  [19]  I.  Lasiecka  and  A.  Manitius  gave  for 
the  first  time  optimal  rates  for  the  AV-scheme.  These  estimates  are  essentially  baised  on 
uniform  (with  respect  to  the  approximation  parameter)  differentiability  of  the  approxi¬ 
mating  semigroups,  which  means  that  the  characterization  of  differentiable  semigroups  in 
[20;  Theorem  4.7]  is  uniformly  valid  for  the  approximating  semigroups. 

In  this  paper  we  develop  a  scheme  using  first  order  splines  which  essentially  has  all  the 
good  properties  of  the  AV-scheme.  In  addition  we  are  able  to  establish  error  estimates 
analogous  to  those  in  [19].  Naturally  uniform  differentiability  is  also  the  essential  basis 
for  our  approach  which  as  far  as  convergence  rates  and  uniform  exponential  stability  are 
concerned  is  motivated  by  the  ideas  of  [19]. 

In  Section  2  we  give  the  basic  ideas  for  the  construction  of  the  approximation  scheme 
using  two  sequences  of  subspaces  leading  to  two  sequences  of  approximating  semigroups 
in  R‘^  X  resp.  dom^I.  We  also  give  matrix  representations  for  the  approximating 
generators  and  the  approximating  input  resp.  output  operators.  Furthermore  we  compute 
the  resolvent  operator  for  the  approximating  generators.  In  Section  3  we  prove  convergence 
of  the  approximating  semigroups  in  both  spaces  by  applying  a  version  of  the  Trotter-Kato- 
Theorem  for  Co-semigroups.  For  the  scheme  in  the  state  space  Z  we  get  also  strong 
convergence  of  the  adjoint  semigroups.  One  should  be  aware  of  the  simplicity  of  the 
consistency  arguments  given  in  Subsection  3.2  for  the  approximating  semigroups  and  their 
adjoints.  Note  that  strong  convergence  also  for  the  adjoint  semigroups  is  needed  for  the 
proof  of  Theorem  7.5.  The  results  of  Section  4  are  fundamental  for  the  rest  of  the  paper. 
In  this  section  we  pro\’e  that  the  approximating  semigroups  are  differentiable  uniformly 
with  respect  to  N.  This  allows  to  use  the  ideas  of  [19]  for  the  proof  of  uniform  exponential 
stability  and  for  proving  rate  estimates.  Uniform  exponential  stability  is  established  in 
Section  5,  whereas  rate  estimates  are  given  in  Section  6.  Implications  of  the  convergence 
results  for  the  approximation  of  the  linear-quadratic  optimal  control  problem  on  the  infinite 
time  interval  are  discussed  in  Section  7.  Finally,  in  Appendices  A  and  B  we  collect  some 
of  the  more  technical  estimates  used  in  Sections  4-6,  respectively. 
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2.  Formulation  of  the  scheme 

In  this  section  we  define  the  spline  approximation  scheme  and  prove  some  basic  proi> 
erties.  Let  =  —kr/N,  k  =  0,...  ,N,  and  =  0,  . 

With  ,  k  =  0,. . .  ,N,  we  denote  the  usual  first  order  basis  splines  on  the  interval  [— r,  0] 
corresponding  to  the  mesh  tg  , . . . 


I  ?(«-<."+■) 
.  0 


for  <  0  <t^' , 
for  <  0  <  t^_j , 
elsewhere. 


Furthermore  we  put 
and 

Eo^'  =  (/„,  0),  =  (0,  t  =  1, ....  A'. 

The  following  spaces  will  be  used  in  the  sequel; 

ir-^'  =  span(f:f'7„,...,F:j;;7„)  C  l7-r,0;R"). 

Z-^'  =  R"  X  W^'  =  span(£-^' . £^')  C  Z. 

A'-"=span(Bo^7„....,Rii:/„)C^f’. 

Z^'  =  C  dom^. 


It  is  convenient  to  introduce  the  ‘'basis  matrices" 


Any  tr  =  (rj.^)  G  can  be  written  as  r  =  =  £•'' col (t/, of - .of  ),  where 

=  col  (of , _ of),  of  6  R",  is  the  coordinate  vector  of  <.5  G  ir'"^  with  respect  to 

the  basis  E^ .  Similarly  any  G  is  given  by  (^  =  B^  b‘’^  with  =  col  (6f  , . . . ,  6f  ), 

G  R". 

The  orthogonal  projections  :  Z  — ►  Z'^  and  are  characterized  in 

Lemma  2.1.  a)  For  {r],p)  ^  Z 


P^7  =  (t^E•"'a•"■), 


where 


Oi'  =  f  / 


V?(s)ds,  fc  =  1, . . . ,  A’. 


b)  For  v  G  -ff  ’ 


Vl-V 


where 


bk 


=  at 


A', 


il-  =  0 . A, 
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Proof:  The  results  follow  by  easy  computations  from  {z  —  r,  )^  =  0,  ^'  =  0, . . . ,  A', 
and  {xj}  —  In)H^  =  0,  A:  =  0, . . . ,  N,  respectively.  | 

Remark.  The  subspaces  were  used  for  the  socedled  scheme  of  averaging  projections 

introduced  in  [2],  whereas  are  the  spline  subspaces  of  [3].  Note  that  P^''  is  not  the 

orthogonal  projection  dom.4  — >  Z^  with  respect  to  the  inner  product  in  Z. 

By  definition  of  the  spaces  Z^ ,  Z^  and  the  operators  A,  B  we  have 

Az^^Z^  fovanyz^eZ^, 

Bi^Z^  for  any  ^eR'". 

For  a  mild  solution  2(t),  t  >  0,  of  (1.3)  as  given  by  (1.4)  we  seek  an  approximation 
iy^(<)  €  f  >  0.  If  2(f)  is  a  strong  (i.e.,  differentiable)  solution  of  (1.3)  then  i(f)  is  in 
general  not  in  dom.4  but  in  the  subspace  generated  by  Az{t)  +  Bu{i).  By  (2.1)  we  have 
Aw^  {t)  +  Bu(t)  G  Z^,  f  >  0.  On  the  other  hand  u'^(f)  is  in  Z^  C  dom>4.  The  above 
consideration  concerning  strong  solutions  of  (1.3)  motivate  to  determine  (f)  such  that 

(2.2)  P^'ii;^’(f)  =  ^P-^'w^\t)  =  Aw^{t)  +  Ba(f ),  f  >  0. 

at 

Remark.  If  instead  of  (2.2)  one  imposes  the  condition 

=  P-^’(.4u;'^'(f)  +  /?n(f )),  f  >  U. 

one  oblcuua  ihe  spline  scheme  of  [5]  which  lacks  a  namber  of  qualitative  properties  one 
would  like  to  have  (see  the  introduction). 


In  order  to  derive  a  differential  equation  for  the  coordinate  vector  of  (f )  we  shall  need 

Lemma  2.2.  a)  P^  restricted  to  Zf'  is  a  bijection  Z^  — »  Z'' .  Its  matrix  representation 
(mth  respect  to  the  basis  of  Z^  and  the  basis  of  Z^  )  is  given  by 


Q^'  = 


/  1  0 
1/2  1/2 
0 


0  \ 


\  0  0  1/2  1/2/ 

b)  A  restricted  to  Z^  is  a  map  Zf''  —*  Z^  with  the  matrix  representation 


(D^  ••• 

•  D^n\ 

(0  ■■■ 

•••  0  \ 

0 

0 

1  -1 

0 

...  0 

N 

+  — 
r 

0 

0 

V  0 

0  / 

VO  ••• 

0 

1  -1/ 

^  j  ^  j^n(  .\ ■+ 1 )  X  ii(A +1 ) 
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where  =  UB^' )  =  Z',.o  («r )  +  5l,A(»)B^ (9)  <16 .  h  =  0, . . . .  A\ 

c)  The  matrix  representation  of  B  considered  as  a  map  into  is 


= 


/^\ 

0 

{o/ 


e  R 


n  ( A'  + 1 )  X  m 


whereas  the  matrix  representation  of  C  restricted  to  is 

C-'^’  =  (C  0  •••  0)  G 

Proof:  For  r  =  (.^(O).y)  G  Zj^  with  ip  —  we  get  (using  Lemma  2.1. a))  z  = 

),  where  i]  =  p{0)  —  b^  and 


,.v 


k  =  1 


.V. 


Part  b)  follows  from  easy  computations  using  the  definition  of  A  and  part  c)  is  trivial.  | 

Let  and  put  z^'(t)  =  P'^'u'^'{t).  i  >  0.  Then  (2.2)  is 

equivalent  to 

(2.3)  z'^'(i)  =  A^z^'(t)  +  Bu{t).  i>0. 

A^  is  an  operator  Z^  — *  Z'"^.  Note  that  (2.3)  is  an  ordinal^’  differentia!  equation  on  Z^ . 
In  view  of  Lemma  2.2  the  matrix  representation  of  is  given  by 

(2.4)  .4 

Let  b^(f)  and  a^{i)  be  the  coordinate  vectors  of  u-^{t)  and  z^(t).  respectively,  i.e.. 
«••'■(/)  =  B^'b^'it)  and  z^'it)  =  E^'a^'it).  Then 

=  Q^'b^\t),  t  >  0, 


and  equations  (2.2)  and  (2.3)  are  equivalent  to 


(2.5) 

Q^'b"'’{t)  =  H^'b^’it)  +  B^'u{t). 

and 

(2.6) 

d^’(i)  =  A^'a^’it)  +  B^'u{t), 

respectively. 


Remarks.  1.  If  A~^  exists,  then  AZ^  =  .  Hence,  the  approximation  scheme  given 

in  this  paper  is  equivalent  to  the  scheme  developed  in  [9]  for  the  approximation  of  the 
spectrum  of  A  using  first  order  splines. 

2.  The  approximation  scheme  developed  here  can  be  regarded  as  a  spline-tau  approxi¬ 
mation  (in  the  sense  of  Lanezos'  so-called  tau-method).  In  [12, 13, 10]  the  Legendre-tau 
approximation  for  delay  systems  was  developed  using  Legendre-polynomials  instead  of 
splines  for  the  definition  of  the  spaces  Z^  and  .  The  Legendre-tau  scheme  was  formu¬ 
lated  in  a  similar  manner  as  above  (see  equations  (2.3)  and  (2.4)  in  [10]). 

3.  Using  different  spline  elements  for  the  definition  of  the  subspaces  Z^  and  one 
can  get  a  whole  family  of  approximation  schemes  using  the  ideas  presented  in  this  paper. 
We  conjecture  that  some  of  them  lead  to  approximations  of  the  characteristic  equation  of 
the  delay  systems  by  functions  which  involve  the  Fade  approximations  of  the  exponential 
function  in  the  main  diagonal  of  the  Fade  table  (see  the  remarks  before  Lemma  2.4) 


For  the  proof  of  convergence  for  the  scheme  presented  in  this  section  we  shall  need  an 
explicit  representation  for  the  resolvent  operators  (\I  —  A^  We  introduce  the  function 


^(t)  = 


1-^/2 
1  -hr/2' 


r  e  C. 


which  is  the  first  non-constant  entry  in  the  main  diagonal  of  the  Fade  table  for  f'’’  (cf. 
[1]).  Next  we  define 

.V 

(2.7)  -r  <  ^  <  0.  A  ^  -2.N77-. 

r=o 


which  is  an  approximation  to  on  [  — r.O]  and  put 

X-^\X)  =  X^-2X/r. 

Proposition  2.3.  aj  .4  complex  number  X  ^  ~2y/r  is  in  a{A^  )  if  and  only  if 

det  A''(A)  =  0. 


b)  Let  X  ^  (?(>!•'■).  A  ^  -2Slr  and  :  =  (r/.E-^a-^')  G  Z^  .  Then 

(A/-^^r’^  =  P'^’(vM0).v). 

where 

rl-(O)  =  A^’(A)->  (rj  +  i:(r^‘)),  =  ej^  UfO)  + 

with 

s  k 

’■>)  =  5^(1  +  »<(rA/;V)) 

A=1  ;=1 

=  ^  E  E I  (»'('■  v.v)*-'  +  ,.(rA/.v)‘-'+')„y 

'  r=i  j=i  ~ 
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Proof:  Given  :  6  we  solve  the  equation  (A7  —  )u’  =  z.  u-  =  G  .  This 

equation  is  equivalent  to  (see  (2.4)) 

(2.8)  {\Q^’  =a^'. 

We  put  =  iQ^  i-e.,  u>  =  with  t/’  =  .  The  definitions  of  Q''  and 

(see  Lemma  2.2)  immediately  show  that  (2.8)  is  equivalent  to 


(2.9) 


N 

Aoq  —  2_^  Die  bie  —  Oq  , 
k=0 

V \  .  Kj  /  t'A 


+  =  +  ^-  =  1 . 


(2.10) 

If  A  ^  —2X/r  then  (2.10)  is  equivalent  to 


b^'  =  p(rA/.V)'^6o^'  +  ^,{1  +  ^,{r\/N))yf>(rX/Xf-^a-^' 

(2.11)  2.\ 

=  cUtk)bo  A-  =  l . A*. 

Therefore  in  case  A  ^  —2X/r  equations  (2.9)  and  (2.10)  are  equivalent  to  (2.11)  and 

N  k 

X^{X)b^  =7  +  ;7y(^  + ///rA/A  ))  ^  ^//(rA/A  )^' 

k=\  j=i 

=  1]  +  I(r-^). 

This  proves  that  A  -2X/r  is  in  a{A^  )  if  and  only  if  det  A-''(A)  =  0.  Moreover,  we 
immediately  see  that  for  A  ^  (7(>l''  )•  A  ^  — 2A/r,  the  element  {XI  —  is  given  as 

stated  in  the  proposition.  | 

Remarks.  1.  Using  the  definition  of  e^  we  see  that 


(2.12) 


(t'^  —  s)E''‘ (s)a'’^  ds 


1. e..  is  the  interpolating  spline  for  the  function  Jg  (^  ~  s)E^  (s)o  ''  ds. 

2.  In  case  A  =  -2X/r  equations  (2.10)  art  equivalent  to 


2#' 

This  together  with  (2.9)  gives 


fA  .-1  M 


^k+l 


7-  =  0,...,A’  -  1. 


,  ^v  -  -«o  “  ■y(^  "  Ty)  i  r  ^ 

"■  ^  A=0 


A'-l 
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Therefore  A  =  —2N/r  is  in  a{A^')  if  and  only  if  det  —  0.  In  case  .4(  )  =  0  and  N 
sufficiently  large  this  is  equivalent  to  det  .4/  =  0,  a  condition  which  is  equivalent  to  the 
existence  of  so-called  small  solutions  for  the  uncontrolled  delay  system  (see  [S]j. 

3.  In  the  single  delay  case  L{^p)  =  Aoip{0)  Ai(p(—r)  we  have 


A^’(A)  =  A/„  -  ,4o  - 


1  -rA/2A^\ 

H-rA/2AV 


N 


‘1 1 


A  #  -2N/r. 


For  the  other  schemes  discussed  in  the  introduction  A^'(A)  in  this  simple  case  is  given 
by  A^(A)  =  XIn  —  Ao  —  f^(X)Ai,  where  /^(A)  is  a  rational  approximation  for 
For  the  scheme  of  averaging  projections  we  have  f^{X)  =  (1  -|-  Xr/N)~^  (see  [2])  for  the 
Legendre-tau  scheme  of  [12]  /^(A)  is  the  A'-th  entry  in  the  main  diagonal  of  the  Fade 
table  for  whereas  for  the  Legendre-scheme  of  [14]  f^’(X)  is  the  proper  rational  entre¬ 
at  position  {N,N  —  1)  in  the  Fade  table  for  For  the  piecewise  linear  scheme  of  [21] 

we  have  /^’(A)  =  (note  that  is  a  Fade  approximant  for  e-"). 

In  case  of  the  scheme  developed  in  [15]  /^(A)  cannot  be  given  explicitely  but  is  obtained 
bv  a  recursion  formula. 


For  later  use  we  state  a  simple  obser\'ation: 

Lemma  2.4.  Let  y  =  ^f=o  l-'Iil  +  J^r^A{B)\dO.  Then  det  A''^(A)  ^  0  for  all  X  £  C  with 
ReX  >  0,  |A|  >  y,  A'  =  1,...  . 

FROOF:  Obviously  |(1  -  rA/2A*)(l  -I-  rA/2A’)“^|  <  1  for  Re  A  >  0.  Therefore  |e3^((?)|  <  1 
for  ReA  >  0  and  -r  <  B  <  0  (see  (2.7)  for  the  definition  of  e'^(B)).  Thus 

( 

J2a,c-^{B,)  + 

1  =  0 

and  the  result  follows  from 


A{B)c^(B)dB 


<  0 


10 


3.  Convergence  of  the  scheme 

In  this  section  we  establish  convergence  of  the  approximation  scheme  presented  in  the 
preceding  section  for  both  state  spaces,  Z  zind  .  Moreover,  we  also  prove  strong  con¬ 
vergence  of  the  adjoint  semigroups  S^{t)*  to  S(f)*  in  the  state  space  Z.  We  first  consider 
convergence  in  Z.  For  that  purpose  a  version  of  the  Trotter-Kato  theorem  as  given  in  [IS] 
will  be  used. 

Let  (1*^,  II  •  11^’)  be  a  sequence  of  Banach  spaces  w'hich  converges  in  the  sense  of  Kato 
to  a  Banach  space  (K,  ||  •  ||)  (see  [18,  Ch.IX,  §3]),  i.e.,  for  each  N  there  exists  a  bounded 
linear  operator  :  1’  — >  1*'^  such  that 

(K  1)  yP^ll  ^  Cl  for  all  N  with  ci  independent  of  N, 

(K  2)  llP'^’yll^  =  llyll  for  all  y  €  r, 

(K  3)  there  exists  a  constant  C2  >  0  such  that  for  all  N  and  all  y  6  1'^  there  exists  an 
X  £  Y  with 

y  =  P^x  and  ||x||  <  C2||y||N. 

Under  these  conditions  the  following  result  is  valid: 

Theorem  3.1.  Let  A  and  be  the  inSnitesimaJ  generators  of  strongly  continuous  semi¬ 
groups  S(i)  on  Y  and  S‘'^(t)  on  1'^,  respectively.  Suppose  that 

(i)  for  some  constants  w-  €  R  and  M  >  1 

II 5 ( t ) II  .V  <  for  all  t  >  0  and  all  N 

and 

(ii)  for  some  A  €  p{A)  fl  0^=1  all  x  ^Y 

lim  ||(A/->l-^’)-*P-"''a--P^’(A/-.4)-’r||.v  =0. 

•N  — 

Then  for  all  x  £  Y 

lim  |15-^'(t)P-’' j  -  P'’'  5(f)x||,v  =  0 
,V  — 3C 

uniformly  for  t  in  bounded  inter^uJs. 

Proof;  The  proof  of  this  theorem  uses  exactly  the  same  arguments  as  those  given  in 
[18,  Ch.IX,  Theorem  2.16]  and  is  based  on  the  equality 

(A/  -  >1'^')-*  (S-’^'COP'^'  -  P^5(0)(A/  -  A)-^ 

-s)((A7-^^')-*P^’  -P^’(A/-^)->)S(s)ds. 

Then  for  any  y  €  dom  A 

(3.2)  lim  |l(A/-^^')-'(5'’'’(f)P'^'-P^5(O)y||.v  =  0 

A  —  oc 


11 


uniformly  for  t  in  bounded  intervals,  where  we  have  used  assumptions  (i)  and  (ii)  together 
with  Lebesgue’s  dominated  convergence  theorem.  Next  we  consider,  for  x  6 


-  A)-^x 

=  (^P"'\XI  -  A)-^  -  {XI  -  A^'r'^P^'^S{t)x 
+  {XI  -  (P^5(0  -  S^{t)P^y 

+  S^(0((A/  -  A^)-^P^  -  P^{XI  -  A)~^y. 

The  second  and  third  term  on  the  right-hand  side  tend  to  zero  uniformly  on  bounded  t- 
intervals  as  N  —*  oc,  because  of  assumptions  (i),  (ii)  and  (3.2).  The  same  is  true  also  for  the 
first  term  on  the  right-hand  side  of  (3.3).  because  for  any  T  >  0  the  set  {S(t)x  1  0  <  t  <  T) 
is  compact.  This  shows  that 


lim  {P^S{i)-  S^{t)P^’)y  =  0 

A'  — oc 


uniformly  on  bounded  t-intervals  for  any  y  G  dom^^l^.  A  density  argument  finishes  the 
proof.  I 

In  our  situation  Z  and  P"^  is  the  orthogonal  projection  Z  — ►  Z  ”^ .  We  have  to 

establish  hypotheses  (K  1)-(K  3)  and  assumptions  (i),  (ii)  of  Theorem  3.1. 

3.1.  Uniform  dissipativity.  Following  an  idea  which  goes  back  to  [2C]  we  shall  use 
equivalent  norms  (|  •  ||a-  on  Z  in  order  to  establish  the  dissipativity  properties  which  imply 
assumption  (i)  of  Theorem  3.1. 

For  j  =  1 . C  -I  define  the  indices  k^'  by  6j  G  _j ).  If  .V  is  sufficiently  large 

(which  we  assume  from  now  on)  then 

1  <  A-f  <  •  •  •  <  A-fli  <  N. 

The  step  functions  {$)  =  g‘^  j)  are  defined  by 


(3.4) 


9k  =  9k^i  +  {orfc=l,...,N-l  with  k  i  {if  |  j  =  1 . ^  -  1), 

9k  =  9k^i  +  ioi  k  =  kj  and  k  =  kj  ~  I,  j  =  I . C  -1. 


It  is  obvious  that 


( 

g{B)  ;=  \im^g^'{d)  =  r  +  ^  ^  X[s^,0)  a.e.  on  [-r.O]. 
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The  inner  products  (•,  ■}x  and  (•,  •)g  are  defined  by 


(('7,</5),(P»V’))i,  =  j  <p{^f4’W9{^)de, 

for  {T),ip),  {p,rl')  G  Z.  Because  of  1  <  g(^),g^{^)  <  ^  +  r,  — r  <  ^  <  0,  the  corresponding 
norms  (( •  ||a',  (( •  ((j  and  ||  •  |(  on  Z  are  equivalent  uniformly  with  respect  to  .V, 

ll.-ll  <  II^IIas  11*11,  <(<’+r)''^||2||  SoTzeZ. 

Therefore  hypotheses  (K  1)  and  {K  3)  are  satisfied  for  (2^.  ||  ■  H.v)  and  (2,  ||  ■  ||j,).  For 
r  =  and  P''  r  =  ft/,  )  we  get 

-  II--BI  = 

This  estimate  shows  that  limy_oc  =  Ikllj  for  ■?  €  Z,  i.e.,  assumption  (K  2)  is 

also  satisfied. 

Lemma  3.2.  Let  t  >0,  be  the  semigroup  onZ*'*'  generated  by  A''^ . 

a)  For  all  z  €  ,  z  — 


')a 


-  4  Hr' I 


0 

L2 


for  X  sufBciently  large,  where  w\\  =  u,-  +  r/X  with 


1  11  .  ^ . 

=  ;^t^n.ax(>lo  +  -4o  )  +  -{(  +  r)  +  -  ^  |*4;1^  +  11-411^3. 

>=1 

b)  For  N  sufficiently  large 

li5'"'(t)|Lv  <  e--'-'  and  jlS^^'CtnU-  <  t  >  0. 

Proof:  Part  b)  is  an  immediate  consequence  of  part  a),  because  {A^'z,z)n  <  wa’||^||^ 
for  all  c  G  Z^ .  In  order  to  prove  a)  let  c  =  Eso^  €  Z^  and  z  =  (0),  )  with 

B  Then 

(3.5)  a^'  =  Q^’b^\ 
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For  elements  zi  =  and  22  =  E‘^  a  simple  computation  yields 

(3.6)  {zuZ2)N=={c^’fR^d‘'^, 

where  =  diag  (l,  , . . . ,  j^g^’)  (8  Using  (3.5),  (3.4)  and  (2.4)  we  get 

(>1^'2,2)a’  = 

Simple  computations  using  (3.4)  and  the  matrices  as  given  in  Lemma  2.2  show 

that 

{A^z,z)^-  =  +  5  E(9l'+,  -  S,'')|f'J'T  - 


k=l 


N 


4- 


it=l 


(3.7)  =  5((  +  r)|6,'-p  - 


fc=l 


*•  J=1  ^  '  '  '  '  / 

+  (C E  (E  )) ‘if’  +  (<f  f  -mi Y.  Bi\«)bA  dO. 

k=l  \j  =  l  ''  d-r  / 

Observing  J3^  (0j)  =  0  for  k  ^  and  k  ^  k^  -  1  we  get 

k=l  ^j=l  ^ 

t 

<  |6;''|  Y  +  |6''’|  1.4, 1  Itjl 

J  =  ] 


J=1 


J  =  1 


Using  Aq^  +  (1  —  \)d^  —  (Aq  +  (1  —  A);d)  >  0  for  0  <  A  <  1  and  q,  /d  >  0  we  get 
N  /  t 


k=l  ^j=\  ' 


(3.8) 


<5ICPEi^^I'  +  5I‘nP 


>=1 

f-1 


;=i 
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Equation  (3.5)  implies 


(3.9)  =  ;(f  E  E  Ed*'*-.!"  + 

k=l  k=l  t=l 

From  -  ^)^k  +  (^  -  ^k-i)^k-i)^  tk  -  ^k-i^  obtain 


llv’''lli>  = 


E(i6?'p + (kHVbt,  +  itfL.i’)  <  ^  E(i‘*-.i'  +  i‘‘'i’)' 


This  implies 


(3.10) 


<  |6^'IMIIl,>llv>'"lll..  <  l|.4|li.|3o'T  +  jlM''lli^ 
<ii-'iiii.i‘?'p  +  5VE(|i>^Mp  +  i3?'n. 


The  estimates  (3.S)  and  (3.10)  together  with  (3.7)  imply 


(A’^x,  =}k  <  -  5^  E  ^  Ed*^'-  d  + 


<wKp--^£(i6f:.p+iir’p) 

*  Ar=l 

+  P  “  4iV  ^  H 

it=i 


Observing  (3.9)  zmd  6^  =  rj  we  get  the  desired  result.  | 
Remarks.  1.  It  is  not  difficult  to  verify  that 


(3.11) 


<  u.'b(0)p  -  -Iblli*  for  all  z  =  (<,5(0),  s?)  €  dom.4. 


2.  In  the  estimate  (3.10)  the  right-hand  side  can  be  replaced  by  (l/2)((.4((^j |i^'P  + 
(r/AN)  Yjk=\ P  +  \^k\^)  which  would  give  the  dissipativity  estimate 

(.4''^  r,  2)a-  <  +  -  ^An,ax(.4o  -l-j4o)-t-f-|-r-|-^  |Ajp  -f 
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This  would  also  be  sufficient  for  b).  On  the  other  hand,  it  was  shown  in  [11]  that  the 
estimate  as  given  in  a)  can  be  used  to  give  a  proof  for  uniform  exponential  stability  of  the 
approximation  scheme  different  from  that  in  the  present  paper  (which  is  based  on  uniform 
differentiability  of  the  approximating  semigroups). 

3.2.  Consistency  of  the  scheme  in  Z.  As  we  shall  see  it  is  considerably  simpler  to 
verify  assumption  (ii)  of  Theorem  3.1  for  and  .4  if  0  €  piA).  In  this  case  obviously 
zilso  0  €  piA^)  for  all  N  because  A^(0)  =  A(0).  In  this  case  we  also  get  assumption  (ii) 
for  (.4^)*  and  A*  without  additional  efforts. 

Lemma  3.3.  Suppose  0  6  p(A).  Then 

lim  \\{A^)-^P^z-P^A-^z\\n  =  0 
N~oo 

and 

lim  =0 

A  —*00 

for  all  z  E  Z. 

Proof:  The  adjoint  of  P‘^  is  the  injection  :  Z^'  — »  Z.  Then  observing  A‘’^(0)  =  A(0) 
and  €^{6)  =  1  we  get  from  Proposition  2.3,b)  that  =  P‘’^  A~^  z  for  c  €  Z‘^ ,  i.e.. 

)~^  —  P^  A~^i''^  and  ((>1'^)*)  ^  =  P'' (.4*  )~^  .  Therefore  for  any  c  €  Z 

WiA-'^'r^P-^’z  -  =  (|P-'^’.4-’(P'’^’--  -  c)|i.v. 

The  result  follows  from  (K  1)  and  P^  z  —*  z.  The  proof  for  the  adjoint  operators  is 
completely  analogous.  | 

The  basic  convergence  result  for  our  scheme  in  the  state  space  Z  is  contained  in 

Theorem  3.4.  Let  S^{t),  t  >  0.  be  the  semigroup  on  Z^  generated  by  i.e.,  S'^{t)  = 
t  >  0.  Then  for  all  z  ^  Z 

lim  S^'{t)P^'z  =  S{t)z  and  lim  S^' {t)* P^' z  =  S(ty z 

A— CSC  iV  — oc 

uniformly  for  i  in  bounded  intervals. 

Proof:  Let  us  define  the  operators  T :  Z  Z  and  :  Z^  — >  Z'^  by 

(3.12)  T(7?,^)  =  (77,0)  {ov{T),^)eZ  and  =  P^’T^a  . 

Replacing  the  operators  A  and  ,4^  by  .4  —  kT  and  .4^  —  kT^  ,  respectively,  just  means 
that  the  matrix  Aq  in  (1.1)  is  replaced  by  .4o  —  ^\In■  Therefore  we  can  choose  k  6  R 
such  that  the  constants  u.’  and  corresponding  to  A  —  kT  and  (where  jz 

and  u!^  are  defined  in  Lemma  3.2),  respectively,  are  negative.  Let  T{t)  and  T^{t)  be  the 
semigroups  generated  by  A  —  kT  and  A^  —  kT^  .  By  Lemmas  3.2  and  3.3  the  assumptions 
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of  Theorem  3.1  are  satisfied  for  T{t)  and  T^\t)  with  A  =  0  in  assumption  (ii).  Hence  for 
all  r  e  Z 


lim  \\T^’{t)P^z  -  P^'T{t)z\\s  =  0 

N  —  oc 


uniformly  on  bounded  f-intervals.  Since  the  norms  ||  •  H;;'  and  ||  ■  \\  are  uniformly  equivalent 
and  ||P'^T’(t)r  —  T(t)r||  — ♦  0  as  .V  — +  oo  uniformly  on  bounded  /-intervals,  we  obtain 


.  (3.13) 


lim  T^(t)P^z  =  T(t)z 

N  —  oc 


uniformly  on  bounded  intervals. 

Returning  to  the  semigroups  5(/)  and  5'^(/)  we  first  observ'e  that  by  the  variation  of 
constants  formula 


(3.14) 


S{f)  =  T{t)  +  K  f  T{t  -  s)TS(s)ds, 

Jo 

+  K  f  T^'{t  -  s)T^S^'is)ds, 
Jo 


for  t  >  0.  Since  T  is  a  finite  rank  operator,  (3.13)  implies 


(3.15)  \\T-^'it)P^’T  -  Tit)T\\ciz,Z)  -  0  as  .V  ^  oc 


uniformly  on  bounded  /-intervals.  From  (3.14)  we  get 

||S(()--  -  <  111(0=  -  T-''(()P-'  =  || 

+  K  f'wnt  -  3)T  -  T'-\t  -  oe'''TlU(z,z,||s«(»)P‘'-’=ll  * 

Jo 

+  K  /'lirif  -0T||^,Z.Z|||S(0.- 

Jo 

By  Gronwall's  lemma  using  Lemma  3.2.  (3.13)  and  (3.15)  we  obtain 

lim  S^\i)P^'z  =  S{i)z 

A  — oc 

uniformly  on  bounded  /-intervals  for  all  z  ^  Z. 

The  proof  for  the  adjoint  semigroups  is  completely  analogous.  | 

Corollary  3.5.  For  z  £  Z  and  u  €  •£'^^^(0,  oo;  R"*)  Jet  z^  (/)  be  the  solution  of  (2.3)  with 
z^(0)  =  P^ z  and  z(t)  be  the  solution  of  (1.3)  with  2(0)  =  z.  Then  for  any  T  >  0 

lim  z^(t)  =  z(t) 

N—oo 


uniformly  for  t  G  [0.  T]  und  uniformly  for  u|(o,t)  bounded  subsets  of  V{0.  T:  R"'). 

Proof:  On  the  basis  of  Theorem  3.4  the  proof  is  quite  analogous  to  the  proof  given  in  [2] 
for  a  similar  result.  | 


17 


Remark.  In  Appendix  B  we  shall  prove  that  (A7  —  A^  )  ^  z  — >  P'' [XI  -  A)  ’2 
as  N  — ♦  oc  for  all  2  6  Z  and  any  A  G  p{A)  (Lemma  B.l).  In  view  of  this  result  the 
introduction  of  the  operators  T  and  in  the  proof  of  Theorem  3.4  is  in  principle  not 
necessary  as  far  as  convergence  of  the  semigroups  S^(t)  in  Z  is  of  concern.  But  we  already 
have  seen  in  this  section  that  the  approach  taken  in  this  section  gives  also  convergence  for 
the  adjoint  semigroups  (<)*  without  additional  efforts.  On  the  other  hand  the  proof  of 
{XI  —  (.4^  )*)  ^ P^’ z  —*  P^{XI  —  A'*)~^z  for  A  ^  0  would  be  very  involved. 

3.3.  Convergence  of  the  scheme  in  We  shadl  prove  convergence  of  the  scheme  in 
using  results  already  obtained  for  the  scheme  in  Z.  For  the  nonhomogeneous  problem 
we  shall  also  need 

Lemma  3.6.  Let  (t),  t  >  0,  be  the  solution  of  (2.2).  Then  there  exists  a  constant  u;  >  0 
(in  fact  Cb  =  (3/2)1|L||£(//»,R’')  +  Ij  such  that  for  all  N 


+3  /  e2-('-)|Bnu(5)l2d5,  t  >  0. 


/■ 


Proof;  Let  lU'^(t)  =  Then 


Since  b^{t)  solves  equation  (2.5),  we  get 


Id 

2di 


1]/- 


=  b^’{tfb‘^\t)  + 


N 


t=l 


k=0  k=\ 


<  (jllillaHMl-)  +  +  ||Bp|u(0P. 


Thus  Gronwall’s  inequality  yields  the  result.  | 


Theorem  3.7.  For  2  =  ((,?(0),<,?)  with  (p  e  and  u  6  Tfj,^(0.  oo;R"')  let  z(t)  be  given 
by  (1.4)  and  w^(t)  be  the  solution  of  (2.2)  with  initial  value  u'^(O)  =  /P/'r-  Then 

lim  llm’^(f)  -  t2(t)llH>  =  0 

A  — OC 
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uniformly  on  bounded  t-intervals. 

Proof:  Since  and  dom>l  supplied  with  the  graph  norm  are  metrically  isomorphic, 
it  sufficies  to  prove  ||>4u’^(t)  —  Az{t)\\z  0  and  (iiD^’(t))(0)  {iz{i)){0)  as  N  —*  oo 

uniformly  on  bounded  t-intervals.  By  definition  of  and  we  have 

(3.16)  A^P^w^{t)  =  AL^P^w^{t)  =  Aw^\i),  t  >  0. 

We  first  consider  the  case  u  =  0.  Using  (3.16)  we  obtain  from  (2.2)  {d/dt)Au:^\t)  = 
A^{d/dt)P^w^{t)  =  A^Aw^it),  i.e., 

(3.17)  ^u>^(t)  =  S^it)Aw^{0),  i  >  0. 

Lemma  2.1  implies 

(3.1S)  P'''.4r(0)  -  .4!i'''’(0)  =  0). 

where  )■  for  some  K  >  0  independent  of  N  we  have  (cf.  [23,  Exercise 

(2.10)1) 

(3.19)  llsr-''  -  9l|t-  <  A-(r/,V)'/"M„, 

for  A’  =  1.2 . This  together  with  (3.17)  and  (3.18)  implies 

||>lu’^’(t)-^--(t)ilz 

<  |i5‘"'(t)l|  |i>lu>‘"''(0)  -  F^U--(0)1U  +  l|S'"’(t)P-^Ur(0)  -  S(t).4--(0)|lz 

<Me-'llI|U(H.,Rn)A'(^)‘/^lblU> 

+  l|S-^'(t)P^’X-(0)  -  S(f)>^;(0)|lz  for  t  >  0. 

By  Theorem  3.4  this  gives 

lim  (/)  -  ^c(t))|z  =  0 

.V  —  3C 

uniformly  for  t  in  bounded  intervals.  Since  r^(t)  =  P‘'^'ui^(t),  f  >  0,  satisfies  (2.3),  by 
Theorem  3.4  we  also  get  lim;^_oo  ||P^’ ui^(t)  —  z{t)\\z  =  0  uniformly  for  t  in  bounded 
intervals.  This  implies 

|(/ui^ (t))(0)  —  (t?(t))(0)|  — >  0  as  N oo 

uniformly  on  bounded  t-intervals.  Note  that  (iw^(t))(0)  and  (iz(t))(0)  are  the  vector 
components  of  P^  [t)  and  z(t),  respectively. 

We  next  consider  the  case  z  =  0  and  u  €  C*(0, r;R'"),  T  >  0.  Equation  (2.2)  together 
with  (3.16)  implies 

p-^  u.-^  (t)  =  /  S^{t  -  s)Bu{s)ds,  t>0. 

Jo 
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Hence  we  have 


Similarly  we  get 


=  A^'P^'w^ii)  =  A^  -  s)Bui3)ds 

Jo 

=  S^’{t)Bu{0)-Bu{t)+  f  S^’{t-s)Bii{s)ds. 

Jo 

Az{t)  =  S{t)BuiQ)  -  Bu{t)  +  [  Sit~3)B{i{s)ds. 

Jo 


Therefore 


A:it)  -  Au'^'it)  =  (5(0  -  5'^'(0)5m(0)  + 


s)  -  S^\t  -  s))Bui3)ds. 


By  Corollary  3.5  this  implies 

hm  \\A=it)-Aw^\t)\\z=0 

A  — oc 

unifornJy  for  t  in  bounded  intervals.  As  in  the  previous  case  (but  now  using  Corollary  3.5) 
we  also  get  lim.v_oc(<ii’^  (t))(0)  =  (ir(<))(0)  uniformly  on  bounded  /-intervals. 

In  case  z  =  0  and  it  G  1^(0,  r;R"’)  w'e  have  the  estimate  (see  Lemma  3.G) 

<  v/36“'|fl|||«||i,.(0.7-;R”).  0<t<T. 

We  have  an  analogous  estimate  for  tz{t).  Density  of  ii^  Z‘(0.r;  R’")  and  a 

simple  application  of  the  triangle  inequality  implies  the  result  for  the  general  case.  | 


4.  Uniform  differentiablity 

In  this  section  we  first  establish  differentiability  of  the  semigroups  5  ''  (t)  uniformly  with 
respect  to  N.  Uniform  differentiability  of  the  approximating  semigroups  S'' {t}  is  funda¬ 
mental  for  our  approach  to  obtain  rates  of  convergence  for  our  scheme  (see  Section  G). 
Also  the  proof  for  the  uniform  exponential  stability  property  of  the  approximating  semi¬ 
groups  (see  Section  5)  is  based  on  uniform  differentiability  of  the  approximating  semigroups 
though  a  different  proof  would  also  be  possible  (see  the  remarks  at  the  end  of  Section  3.1 ). 
For  a  fixed  real  number  b  >  r  and  numbers  a  >  0,  u;  >  0  we  define 

Sa,w  =  {A  G  C  1  llmAl  >  and  ReA<u;} 

and  the  “exponential  sector" 

=  { A  G  C  I  Re  A  <  u.’  and  A  ^  }  • 

Theorem  4.1.  a)  Let  u:  >  h  in  Lemma  2.4),  b  >  r,  0  <  d  <  1  and  put  oi  = 

u>  -i-  jin  Y^.  Then  for  all  a  >  max(ooi )  (oo  as  defined  in  Lemma  A. 5) 

(T{A^')CSa.^  for  ally 

and  there  exists  a  positive  constant  ci  such  that 

ll(A/- <ci|ImA| 


for  all  A  €  Ta,.'  and  all  .V. 

b)  For  any  compact  subset  K  of  p{A)  there  exist  positive  constants  C2  and  .Vo  such  that 
K  C  p(A-^')  for  N  >  N2  and 

IK  A/  -  A^'  )■’  II  <  02  for  A  G  K  and  .V  >  .Vj . 

Proof;  By  definition  of  oi  and  A*o  we  have 

(4.1)  A'oe““*  <  - -  and  |ImA|  >  - - -  for  A  e  Ua.^-. 

We  first  prove  that  A  G  p(A'^  )  for  all  A  €  Uo.w  and  all  A'.  According  to  Theorem  2. 3, a) 
this  is  equivalent  to  det  A''(A)  ^  0  for  A  6  Ua,w  Lat  a  >  max(ao,ai).  Then  Lemma  A. 5 
and  (4.1)  imply 

lA(e^^7„)|  <  7^max^jMAr/A7l'‘'  <  w7voe""^|Im  Aj 
<(1-<^)|A(  for  A€E„,^,  A^=  1,2,...  . 

Therefore  |A^(A)|  >  ^|A|  for  A  €  So,.-,  N  =  1,2,...,  and  (taking  the  Xeuman  series  for 

A'''(A)-') 

(4.2)  |A^'(A)-M  <  for  A  G  .V  =  1.2 . 

d|Aj 
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This  proves  a).  Note  that  by  Lemma  2.4  any  A  with  ReA  >  ->  is  in  )  for  all  N. 
According  to  Theorem  2.3,b) 


{XI  -  =  P^’(V’(0),  v'>)  for  all  r  6 


where  ?/’  is  given  in  Theorem  2.3  and  2  =  a  .  Using  Lemma  A. 5  we  get 

(/i  =  fi{rX/N)) 


j=i 
k-l 


2N 

1  r 


(4.3) 


i=i 

k 


J=1 


^/U-S. 


J=1 


and 

(4.4)  |I(^-"’)l<-r*/'(l-<5)||^|U5|ImAl 

u.' 

for  all  A  e  Eq.w  and  all  .V.  The  estimates  (4.1).  (4.2)  and  (4,4)  imply 

1W0)|  <  ^(|„„''|+ir'/"i^|k-|U,|lmA|) 

(4.5) 

C|A|  \ 


^|A| 

.  1-6/1 


for  all  A  £  Ea,w<  2  £  and  jV  =  1,2, . . .  .  Using  (4.1)  again  we  also  have 

(1-6)2/! 


(4.6) 


lv-(0)|<^-^(-  +  r>/2)||c||  |ImA| 


for  all  A  £  z  £  and  A  =1.2 . From  (4.3)-(4.5)  and  Lemma  .A. 5  we  obtain 

(4.7)  <  l^(Ljl{l  +  r>/2)  +  r>/2)||,||  |in,  A|.  k  =  \ . .V. 

90 


for  all  A  e  r  e  and  -V  =  1.2,...  .  Let  cj  =  ^  +  r^/^)  +  r^/^) .  Then 

(compare  also  (3.9))  by  (4.6)  and  (4.7) 

l|j’''(v(0i.e)|p  <  |V'(0)|'+  ^ 

it=l 

<c?(l+r)||z(j2|ImA|2 
for  all  A  E  Lla.w,  2  G  and  A’  =  1,2, ...  . 

In  order  to  prove  b)  we  use  Corollary  A. 4  to  see  that  K  C  p{A^  )  for  N  sufficiently  large, 
say  N  >  N2  (compare  Proposition  2.3,a)).  We  also  see  that  A'^(A)“^  -+  A(A)“^  uniformly 
on  A',  which  implies  |A^(A)“^|  <const.  for  A  €  K  and  N  >  N2-  Furthermore  Corollary 
A. 3  shows  that  sup_,.<g<Q  |e^(^)|  =  maxo<j<A'  |^(rA/A'’)|-'  <  const,  for  A  6  A' and  aJl  N. 
Therefore  (4.3)  reduces  to  |7"^(t^  )j  <  const.|j<p|]£,2  for  A  €  A'and  all  X.  The  rest  of  the 
proof  uses  the  corre.'^pondingly  simplified  versions  of  (4.4)-(4.7)  in  order  to  get 

|!P'^  (!.'(0).  t')||‘  <  const. ||r||^  for  A  G  K  and  N  >  A'l .  | 

It  is  easy  to  see  that  for  arbitrt  n'  o.u,’  G  R  and  b  >  r  we  have 

(4.8)  <  A'oe~“'’lImAl,  -r  <  ^  <  0,  A  G 

If  we  choose  u;  >  a  >  max(ao-0))>  then  analogous  computations  as  in  the  proof  of 
Theorem  4.1  show  that 

||(A/ ->!)■' II  <  cilImAI  for  all  A  G 
Let  r=riUr2Ur3.  where 

r,  =  {A  G  c  1  ReA  <,.>  and  Im  A  = 

=  {A  G  C  I  ReA  =  and  lImA|  <  et"— '>*}. 

Fa  =  {A  G  C  I  A  gFi). 

The  following  theorem  i.*;  an  immediate  consequence  from  [20;  Thm.  4.7]. 

Theorem  4.2.  a)  For  all  A’  and  all  b  >  r 


(4.10) 

b)  For  all  b  >  r 

(4.11) 


(4.12) 


S^(/)=;^  /  d\,  t>2k 

Jr 

A^'s^'(i)=^  f  d\.  t>: 

2-t  Jr 


>  36. 


c*'(A7- ^)-’dA,  />26, 


Xc’^^iXI  -  A)-^dX.  t>3b. 


Reniark.  The  contents  of  Theorems  4.1  and  4.2  explain  why  we  labeled  our  scheme 
“uniformlv  differentiable". 
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5.  Uniform  exponential  stability' 

In  order  to  establish  the  property  of  uniform  exponential  stability  for  our  scheme  we 
shall  use  convergence  of  the  resolvent  operators  {XI  —  .  Our  approacli  is  the  same 

as  in  [19]. 

Theorem  5.1.  Let  Uo  =  sup{J?eA  |  A  6  cr{A)}.  Then  for  any  e  >  0  there  exists  a 
consteMt  M  =  M(e)  such  that 

||S^'(t)ll  <  t  >  0, 


for  N  sufRciently  large. 

Proof;  By  assumption  we  have  A  G  p(A)  for  ReA  >  uq.  Corollary  A. 4  together  with 
Proposition  2. 3, a)  implies  that  for  N  sufRciently  large  A  G  p{A^)  for  ReA  >  cjq  +  c- 
Therefore  we  can  choose  a;  =  wo  +  e  in  the  definition  of  the  path  P  for  the  representation 
(4.9).  This  is  a  consequence  of  Cauchy’s  theorem,  because  {XI  —  -4)~^  and  {XI  —  A^)~^ 
for  N  sufficiently  large  are  analytic  in  ReA  >  u,’o  +  e.  Furthermore  we  choose  b  =  5r/4. 
From  (4.9)  we  get  the  estimate 

I1S'U)||  <  /  £''>'"||(A/-.4’'')-M|WA| 

^  /  ||(A/->f^’)-M||dA|  =:  +  J2  fort>2?>. 

Jr^ 

Using  Theorem  4.1. a)  the  estimate  for  Jj  is 

Ji<-f  c'’^^^|ImAi|dA|.  /  >  2b. 

^  Jri 

A  parameter  representation  for  Fj  is 

(5.1)  A(r)  =  r -f- -oc  <  r  <  +  e, 

with  |A(t)|  =  (1  +  Therefore 


J]  <  const 


/-’O  +  f 


.T(l-2b) 


dr  < 


const . 
t~2b' 


(-•o+Oi  <  Cf(-.o  +  0(  fyj.  f  >  3,._ 


where  c  is  appropriately  chosen  (note  that  2b  =  5r/2  <  3r). 

The  path  Fj  is  a  compact  subset  of  p{A).  By  Theorem  4.1.b)  we  have  for  A  G  F2  and 
N  sufficiently  large 

J2  <  const. for  t  >  3r. 


which  together  with  the  estimate  for  Jj  shows  that 

(5.2)  ||S‘’^(t)||  £  const. for  i  >  3r. 

By  Lemma  3.2  (||5'‘'‘ (t)|l)yy_j  ^  is  uniformly  bounded  on  [0,3r].  This  together  with 

(5.2)  prov'es  the  result.  | 

Remark.  If  u.’o  <  0  then  Theorem  5.3  states  that  the  approximating  semigroups  S^’{t) 
are  exponentially  stable  for  -V  sufficiently  large  with  a  uniform  decay  rate,  which  can  be 
chosen  arbitrarily  close  (from  above)  to  the  decay  rate  of  the  semigrouj)  S(i). 
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6.  Convergence  rates 


In  this  section  we  establish  rates  of  convergence  for  our  scheme.  The  basic  idea  for  the 
approach,  namely  to  use  uniform  differentiability  of  the  approximating  semigroups,  was 
introduced  by  I.  L^lsiecka  and  A.  Manitius  in  [19]  for  the  scheme  of  averaging  approxima¬ 
tions.  We  shall  prove  two  kinds  of  rate  estimates:  a)  Optimal  rates  on  t  >  0  for  sufficiently 
smooth  initial  data,  b)  Optimal  rates  on  i  >  to  >  0  for  general  non-smooth  data. 

For  fixed  w  >  7  (7  eis  defined  in  Lemma  2.4)  and  e  >  0  we  choose  b  £  {r,r  +  e/7)  and  a  = 
m2Lx(ao,ai)  (oo  defined  in  Lemma  A. 4  and  ci  in  Theorem  4.1, a))  in  the  definition  of  Sa,w 
Let  Wo  =  sup{ReA  |  A  G  For  any  u>i  £  (wo,w]  the  path  Fu,,  =  Fi.u;,  U  Fj.u-i  U 

is  defined  by 


Fi.w,  =  {AgC  lReA<wi  and  Im A  = 

F2,....  =  {A  €  C  I  ReA  =  wi  and  [ImAl  < 

Fa,.,  ={AgC|AgFi.,..}. 

The  following  assumption  will  be  used  in  some  of  the  estimates: 

(H)  The  delays  Oj,  j  =  are  commensurate  (i.e.,  for  f  >  0  and  natural  numbers 

kj  we  have  9j  =  A'jf.  j  =  1 _ ,()  and  are  contained  in  the  set  of  meshpoints  , 

A'  =  1, . . . ,  A*  (which  is  the  case  for  N  =  mk(,  m  =  1.2,...). 

The  first  result  on  convergence  rates  is  concerned  with  norm  convergence  of  the  approx¬ 
imating  semigroups  (with  rates  0(1/A’)  resp.  0(1/A’’/^))  for  f  sufficiently  large  and  with 
convergence  uniformly  on  bounded  intervals  if  initial  data  are  in  dom>A^. 

Theorem  6.1.  Let  a  =  1  if  (H)  is  satisfied,  otherwise  a  =  1/2. 

a)  For  any  e  >  0  there  exist  positive  constants  Cj,  C2  and  A'l  such  that  for  A’  >  A’j 


(6.1) 

f  >  5r  -1-  e. 

and 

(6.2) 

llI^'S‘''(t)P' 

,  f  >  6r  -F  e. 

b)  For  any  T  >  0  there  exist  positive  constants  Cj  and  N2  such  that  for  N  > 

(6.3) 

||S^(t)P-'^ 

■z-P^'S(t)z|l<C3(^)°llzlUn>.4’, 

0  <  t  <  r. 

for  all  z  £  domA'^  ('ll^lldom =  \\  =  \\  +  M^-ll  for  z  £  dom^A^). 

Proof:  For  c  >  0  and  sufficiently  small  let  0,7  =  wq  +  e.  By  Theorem  4.1, a)  we  have 
)  C  5a,w  for  all  A^  From  Theorem  5.1  we  infer  that  ReA  <  wi  for  all  A  G  a{A^') 
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and  N  sufficiently  large,  say  N  >  .  By  Cauchy’s  theorem  we  can  replace  F  in  Theorem 

4.2  by  Fu,!  for  N  >  N^.  Using  (4.9)  and  (4.11)  we  obtain 

||P'''5(0-5-''’(t)P^’||  <  ^11  /  e^‘(p'^(A/->l)-'  -(AJ->l^')-*P'"’)dA|| 

3,4*/ I 

+  ^11  /  c^'(p^(AJ ->!)-'  -(XI  -A^’r^P^’]  dAll 

•'t'  J.wi 

Using  the  parameter  representation  (5.1)  for  Fi^i^,  and  Lemma  B.l,a)  we  get  (note  that 

56  <  5r  +  e) 

AiMjjy  f  e""‘|ImAnrfA| 

1  ,4*/  I 

r  ^  « 

<  const. (  —  )  j  <  const. (  — 

for  t  >  5r  +  e  and  N  >  N^.  Using  Lemma  B.l,b)  we  obtain  analogously 

J2  <  const. 

for  t  >  5r  +  €  and  N  >  Ni.  Combining  the  estimates  for  J\  and  J2  we  get  (C.l). 

For  the  proof  of  (6.2)  we  first  note  that  S^’{t)P^z  -  S{t)z  =  {i{L^  S''^' {t)P^  z  - 
S(t)s)(0), . . . )  for  r  €  dom,4.  Therefore  (6.1)  implies 

(6.4)  |,(I'"’S''’(f)P-'’r  -  S(f)--)(0)|  <  c,(^)%<-'>+'>' 

i  ' 

for  t  >  5r  +  c  and  N  >  A*i .  It  remains  to  estimate 

-  Al5(0IU(^.Z) 

(«-5)  <  -  P^'AS(t)Hc,z.z> 

+  ||P'"a(5(()  -  AS(t)Uz,zy 

From  (4.10)  and  (4.12)  we  immediately  get 
S^' (t)P^’ -  P^' AS(t) 

=  7^-  J  -  A^'r'^P^'  -  P^'iXI  -  A)-^^  dX 

for  t  >  3b  and  N  >  Xi.  Estimates  analogous  to  those  leading  to  (6.1)  (using  Lemma  B.l) 
give 

(6.6)  \\A^S^’it)P^  -  P^U5(/)|U(^.^)  <  const.(^)%<-o+'>' 
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for  t  >  6r  +  e  and  N  >  N\.  Note  that  in  the  estimates  we  pick  up  an  additional  fac¬ 
tor  |A|  which  is  bounded  by  “const. |Im A|”  on  We  still  have  to  estimate  ||(P‘^  — 

I)AS{t)\\c(Z.Z)-  Let  t  >  2b.  Since  S{t)  is  differentiable  for  t  >  r,  we  have  >^S(2i)r  €  dom^ 
and  therefore 

\\P^’AS{t)z  -  ASit)z\\  =  llP^’5(t  -  2b)AS{2b)z  -  5(t  -  2b)AS{2b)z\\ 

^  <  const. ^llt5(t  —  26)^5(26)2:llH»- 


Here  we  have  used  that 

||P^'((^(0),  c^)  -  (<p(0),<p)||  <  const. for  <p  £  H\ 

(compare  [19;Proposition  A. 3]).  By  differentiability  of  S{t)  we  get  from  [20;  Lemma  4.2] 
that  AS(2b)  and  A'^S(2b)  are  bounded  operators.  Thus 

||/S(t  -  2b)AS(2b)z\\L.  <  |lc|l 

and 

||t^S(t  -26MS(26)c|li..  =  ||/5(t  -26M''S(26)zlU,  <  e<“'‘>+‘><'-2*)l|>l2S(26)||  |lrl| 
i.e..  ||tS(t  —  26)^5(2?))- IIh'  ^  const. which  together  with  (6.7)  gives 
(6.8)  \\P^"''AS(t)z  -  ^^(O^li  <  const. 

for  all  t  >  2b  and  N  >  N\.  The  estimate  (6.2)  now  follows  from  (6.4)-(6.6)  and  (6.8). 

In  order  to  prove  part  b)  of  the  theorem  we  choose  A  6  p{A).  Then  Lemma  B.l,b). 
Lemma  3.2.b)  and  (3.1)  imply 


||(A/  -  A-^'r^  {S^\t)P^'  -  P'’''5(0)^'||  <  const.(^)“||--||dom.4 

for  all  t  £  [0.  Tj.  c  £  dom^4  and  .V  sufficiently  large.  This  together  with  (3.3).  Lemma 
B.l,b)  and  Lemma  3.2,b)  implies 

||(P'''5(0-5''(()e''')(A/-.4r'.-|| 

<  c(^)°||5(()s||  +  const.(^)°||s||d„„.«  +  const. 
for  0  <  t  <  T,  N  sufficiently  large  and  r  €  dom^.  The  last  inequality  shows 

||/>■'■5(^)c  -  <  const.(^)’“||c|Um.4. 


for  0  <  t  <  T.  A’  sufficiently  large  and  z  £  dom^^.  | 

The  rest  of  this  section  is  devoted  to  optimal  convergence  rates,  i.e.,  rate  estimates  of 
the  form  (9(1/A'^)  and  (9(1/A’^'^^).  respectively. 
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Theorem  6.2.  Let  (3  =  2  if  (H)  is  satisfied,  otherwise  j3  =  3/2. 

a)  For  any  e  >  0  there  exist  positive  constants  c\,  C2  and  A'l  such  that  for  X  >  N\ 

(6.9)  \\S^’it)P^z  -  P^’S{t)z\\  <  t>er  +  e. 

and 

(6.10)  -  Pi^iSCO^IIh.  <  (  >  7r  +  e, 

if  z  6  dom  A. 

b)  For  any  e  >  0  there  exist  positive  constants  C3,  C4  and  N2  such  that  for  N  >  N2 

(6.11)  IIS''’)/)/’-".-  -  P'''S(()--||  <  C3(^)%'“"+'>'(k|  +  ||i2||„-..-).  1  >  6r  +  j. 

and 

(6.12)  ||,I-’''S''’(()P'''--  -  P,'''i5(().-|Ih.  <  C4(;^)"e'““+'>'(|>;|  +  ),  (  >  7r  +  f, 

for  all  z  =  with  r,0: R"). 

c)  For  any  T  >  0  there  exsist  positive  constants  C5  andNs  such  that 

(6.13)  ||S-’''(()P''--  -  P'''5(().-||  <  C5{^)*||2|U„,  J.,  0  <  (  <  r. 

for  all  z  €  dom  .4^  and  N  >  A'3.  For  z  €  dom^l'*  we  can  replace  the  right-hand  side  of 
this  estimate  by  C5(^)  l|rl|dom 

Proof:  As  in  the  proof  of  Theorem  6.1  we  get  for  A’  sufficiently  large,  A*  >  A'2. 
ljP-’^'S(t)z  -  S'^'{t)P-^'zjl  <  ;^|1  /  e^'(p-^’(A/->l)-^--(A/->;"')-'P-^  c)dA|l 

+  ^11/  e'''(p-"(A7  -  Ar'z  -  {XI  -  A' r'P^'c)  dA|| 

=:  ^(P.  +■^2). 

Using  the  parameter  representation  (5.1)  for  and  Lemma  B.2,a)  we  get  (observe  also 
66  <  6r  +  e) 

Ji  <  2c(^)^||rl|doni^  /  e'^'^IImA^dAI 

^  ’'Flu,, 

<  const.(-^)‘^||r||dom.4  /  c’’''"®*’^dr  <  const. (^)^e*‘*'®-*-'^'||rl|dom.4 

J  —  oc  *' 


for  <  >  6r  +  e  and  N  >  N2-  Using  Lemma  B.2,b)  we  obtain  analogously 

J2  <  const.(^)^e<‘^®'*'*^'||2|jdom-4 

for  t  >  6r  +  e  and  N  >  N2.  If  ^  r,0; R")  then  in  the  estimates  for  J]  and  J2 

we  can  take  ^  =  2  and  1?7|  +  ||v5||vv>.~  instead  of  HzHdom^-  This  finishes  the  proof  for  (6.9) 
and  (6.11). 

Since  [iL^' z  -  Pl^ iS{t)z)(Q)  equals  the  vector  component  of  S^{t)P^z  - 
P^' S{t)z,  we  get  from  (6.9)  the  estimate 

-  Pi^i5(<)^)(0)|  <  Ci(^)^(‘"‘>  +  ‘>'|U-||dom^ 


(resp.  <  ci(-^)^e^‘^'>‘'''*'(|//|  +  |l^||ir>.«))Ior  f  >  6r4-e  and  A'' >  N2.  It  remains  to  estimate 
\\iAL^’ S^\i)P^' z  -  iAPl^S{i)z\\i^.  <  \\A^S^’it)P^z  -  APi^' S{f)z\\. 


where  we  have  put  P^^  =  i  i.e.,  P/'’  is  the  orthogonal  projection  dom^  — >  Zf  . 

Since  A^  and  are  bounded  operators,  we  get  from  (4.9)  and  (4.11) 

A^'S^\i)P^'z  -  APi''Sit)z 

e^^(A^\\I-A^')-^P^'z-APl'\\r-A)-^z^d\ 

for  t  >  2b.  Easy  computations  using  also  {A'^  P‘''  =  P^  A~^  P'"^  and  P'^' ✓4Pj'''  = 

Pj^'  yield 

A^’(\i  -  A'^'r^p^'  z  -  APi'ixi  -  Ayh 

(6.15)  =  A^’iXI  -  A^'y'^  {AP,^'  -  P^'A){\I  -  Ayy 

+  xA^’ixi  -  A^’y'^p^’ii  -  p!'’){xi  -  Ay^z. 


From  A^{XI  —  >4^ )“’  =  X{XI  —  A’^'y^  —  I  and  Theorem  4.1, a)  we  get 

(6.16)  (A/ ->!•'' )~*P'^  11  <  const. [ImAj^  for  A  G  Sa,u>  and  A"  >  A'2. 
With  the  notation  used  in  the  proof  of  Lemma  B.2  we  have 

(6.17)  (/  -  P,^’)(A/  -  Ayy  =  (V’(0)  -  -  r/>^)  =  (0,  V’  -  V^’)- 

An  easy  computation  (see  Lemma  2.1  for  the  definition  of  P^  and  P^^  )  shows 

{AP,^'  -  P^’A)(il'{0).xl’)  =  -  v''),0). 
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This  and  (6.15)-(6.17)  imply 

||^''^(A7  -  A^'y'^P^'z  -  APl^iXI  -  Ay^zW 

<  const. |A1  jim  A]^  ’’I'  —  4’^  ) 

<  const.|ImAp(lI(T/)^'  -  V>)|  +  ||v'’  -  lUO 


for  A  €  Ea.u;  and  N  >  N2-  The  estimates  used  in  the  proof  of  Lemma  B.2  show  that 

11.4^(A/  -  A^yP^'z  -  APl^iXI  -  Ayz\\ 

<  const.jlm  A|®(^)^l|2lldon,.4  for  A  e  Sa.w  and  A"  >  .¥2- 

If  r,  0;R”)  then  0  =  2  and  jl-jjldom^t  can  be  replaced  by  It;]  +  Hi^Hivi.oo .  An 

emalogous  proof  shows  that 

WA^’ixi  -  A^'yp^'z  -  Apyxi  -  Ayz\\ 

^  ^  <  const. (^)^||2l|dom.4  for  A  G  r2,^.i  and  N  >  A'2- 

Estimates  (6.18)  and  (6.19)  together  with  (6.15)  yield 

||_4A’5N(,)p.v.  _  ^pA5(().||  <  const.(^)"e<--"+»'||rll„<.„^ 


for  t  >  7r  +  e  and  N  >  A’2  with  the  obvious  modification  in  case  ^  G  Tr''^(— r, 0;  R"). 

For  the  proof  of  part  c)  we  choose  A  G  p{A).  Then  Lemma  B.2,b),  Lemma  3.2.b)  and 
(3.1)  imply 


1|(A/  -  Al-'r'  (S'(()P-'’2  -  P-''S(()--)||  <  const.(.(:)‘*||.-||j.„,^= 

for  0  <  t  <  r,  2  G  domv4^  and  N  sufficiently  large.  This  estimate  together  with  (3.3)  and 
again  Lemma  B.2,b),  Lemma  3.2,b)  implies  (6.13).  The  result  for  r  G  dom  is  clear  if 
we  observe  that  |7?|  +  ||.^||ivi.c^  <  const.|lz||doni.4*-  I 

Remarks.  1.  The  result  of  Theorem  6.2,c)  shows  that  the  vector  component  of  5'^  {t)P‘’^  z 
converges  to  the  vector  component  of  S{t)z  (which  is  the  solution  x{t)  of  (1.1))  with 
rate  1/N^  uniformly  on  bounded  t-intervals  for  sufficiently  smooth  initial  data  (e.g.,  z  G 
dom>t^).  Note  that  the  vector  components  of  S{t)z  and  P^  S{t)z  coincide. 

2.  The  condition  z  G  dom.A  for  the  estimates  (6.9)  and  (6.10)  can  be  replaced  by  the 
condition  z  =  (77,1^),  9  G  TT^'^(— r,0; R"),  as  can  be  seen  from  (B.14). 

3.  Since  S{i)  is  compact  for  t  >  r,  the  estimate  (6.1)  also  implies 

-  5(()|U|z.2|  ^0  as  JV  ^  00 


for  t  >  br  +  €. 

4.  One  should  observe  that  the  estimates  (6.10)  and  (6.12)  are  estimates  for  the  Tf ’-norm, 
i.e.,  involve  the  derivative  of  iL^  {i)P''^  z  —  Py  S(t)z. 
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7.  The  linear  quadratic  optimal  control  problem 
We  consider  the  following  control  problem: 

For  given  initial  data  z  =  E  Z  minimize  the  cost  functional 

(7.1)  J{u,z)=  (ly(t)p  +  lu(#)|^)  dt 

Jo 

over  u  €  £^(0,  oo;R'")  subject  to  (1.1)  or,  equivalently,  to  (3.1). 

Definition  7.1.  a)  The  pair  (A,  B)  is  stabilizable  if  there  exists  an  operator  fC  €  jC(Z,  R”*  ) 
such  that  A  —  BK  generates  an  exponentiai/y  stable  semigroup  on  Z. 

b)  The  pair  (C,  A)  is  detectable  if  there  exists  an  operator  Q  6  £(R^,  Z)  such  that  A  —  QC 
generates  an  exponentially  stable  semigroup  on  Z. 

The  following  theorem  is  well-known  (see  [24]). 

Theorem  7.2.  Let  (A.B)  be  stabilizable  and  (C,A)  be  detectable.  Then  the  algebraic 
Riccati  equation 

(7.2)  (A*n  +  nA-nBB*n  +  rC)z  =  0,  zedomA, 

has  a  unique  self-adjoint  and  non-negative  solution  IT.  Moreover,  the  operator  A  —  BB*Il 
generates  an  exponentially  stable  semigroup  Tit),  i  >  0,  on  Z.  The  optimal  control  for 
(7.1)  is  given  by 

t5(t)  = -rnr(t)z,  t>o. 

Remark.  Without  loss  of  generality  we  can  assume  in  Definition  7.1  that  C  dom  A* 

and  ^R^  C  dom  .4.  In  fact,  the  solution  11  of  (7.2)  satisfies  IIZ  C  dom  .4*  (see  [7]).  so  that 
nSR'”  C  dom  .4*.  Similaily  the  dual  Riccati  equation 

(AE-h^LA^ -ECXi:-hBB*)z  =  0,  zedomA*, 

has  a  unique,  self-adjoint  and  non-negati'''e  solution  E  such  that  SZ  C  dom>l  and  A—EC*C 
generates  an  exponentially  stable  semigroup  on  Z  (cf.  [22]).  Thus  EC*R^  C  dom .4. 

Let  B^  =  B  and  =  C\zn.  The  following  result  is  concerned  with  stabilizability 
of  (A^ ,B^ )  and  detectability  of  (C^ ,A^),  both  uniformly  with  respect  to  N . 

Theorem  7.3.  a)  Suppose  that  (A,B).  Then  there  exist  constants  Mi  >  1,  u.’]  >  0  and 
a  sequence  AC^  G  £(Z^,R'")  such  that  ||^^1|  is  uniformly  bounded  and 

^  f>0, 

for  N  sufficiently  large. 

h)  Suppose  that  (C,A)  is  detectable.  Then  there  exist  constants  M2  >  T  <->-'2  >  0  a 
sequence  0^  €  £(R^,  Z  '^ )  such  that  ||  is  uniformly  bounded  and 

t>0. 
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for  iV  sufficiently  large. 

Proof;  a)  Let  (A,S)  be  stabilizable.  Then  there  exists  an  operator  /C  G  £(Z,R’”) 
such  that  A  —  BfC  generates  an  exponentially  stable  semigroup.  We  put  K"''  = 

=  1, 2, . . .  .  Then,  for  z  G  , 

-  B^K^)z  =  A^z  -  BKz  =  AL^z  -  BKP^’L^'z 
=  {A-BfCP^)L^'z, 

i.e.,  the  semigroup  ^  is  the  approximatiilg  semigroup  for  the  delay  system 


x(t)  =  LN{Xt), 

which  corresponds  to  ^  —  BKP^ .  Let  /C(r/,9)  =  A'o??  +  f^^K(0).p(0)  dO  for  (r],ip)  G  Z, 
where  Kq  G  R"’^"  and  K(-)  G  L^(— r,0;R"’^").  An  easy  computation  shows  that 


K^P'"^  (rf,  lyp) 


K^' (0)^(0)  dO. 


where  =  '£j=i  A'f  with  K^'  =  {N/r)  jl^''K{9)df>.  Therefore 

/•o 

Ia-(^)  =  L{^) - J5A'o9(o) -b J  K^\e)^{e) do. 

The  delay  system  corresponding  to  >1  —  BIC  is  x(t)  =  L{xi).  where 

L{^}  =  L{.p)-BKo^{0}-B  J  KiOM9)dB. 


We  choose  7  as  in  Lemma  2.4  corresponding  to  L  and  u.’  >  7,  6  >  r,  0  <  6  <  1  as  in 
Section  4.  The  estimate 

|iN(eJ7„)|  <  |£(eJ7„)|  +  |iN(e!:'/„)  - 

<  (1  -  i)|A||B|  jjKW  -  K'-'m  |£a  (6)1  de 

<  (1  -  «)|A|  +  |B|r>'^|lA'  -  A^'lli.  max  iK(Ar/.V)l* 

0  <  it  <  A 

<  (1  -  6)\X\  +  |R|r*/2|lA'  _  A'^IU, A'oe-“‘|Im  A| 
for  A  G  shows  that  for  6  G  (0,6)  there  exists  an  Nq  such  that 

I^A’(eA'/n)|  <  (1  -  6)|A|  for  A  e  Ea.u.  and  N  >  Aq. 
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Using  this  estimate  we  can  establish  Theorems  4.1  and  4.2  for  the  operators  .4^  —  , 

N  >  Nq.  But  then  also  Theorem  5.1  is  valid  for  A  —  BK,  and  —  B^ ,  N  >  No,  which 
proves  the  result  on  stabilizability. 

b)  Let  {C,A)  be  detectable.  Then  there  exists  an  operator  Q  €  C{RP,Z)  with  C 
dom>l  such  that  A  —  QC  generates  an  exponentially  stable  semigroup  T{t).  t  >  0,  on  Z. 
Since  T(t)  is  not  the  solution  semigroup  of  a  delay  system  the  proof  in  this  case  will  be 
somewhat  different  compared  to  that  for  a). 

For  iV  =  1,2, . . .  we  define 

e  /:(R^z^). 

Since  ffp  E  dom>l  for  all  y  G  R^,  there  exists  a  function  G  €  H^(—r,0;  R"^^)  such  that 

(7.3)  Gy  =  (G(O)y.Gy)  for  all  y  G  R^. 


Using  the  definition  of  and  Lemma  2.1  we  get 

(7.4)  G^'y  =  (G(0)y,G^y)  for  all  y  G  R^ 


where  G^ 
given  by 


=  (1/2)  Pj^  (G(t^)  +  The  matrix  representation  G^  of  G'^  is 

=col(Go^'....,G’i^;)  G 


where  Gq  =  G(0)  and  G^  =  (l/2)(G(t^_j )  +  G(t^' )),  k  =  I, . . . ,  N .  From  (7.3)  and  (7.4) 
it  is  easily  seen  that 


A’ 

lie''  -ffl!’  <  -  G(«)|  +  IGdi'L.)  -  G(«)|)  <  ■•*(^:G) 


where  a.’(f;  G)  =  siip|^_,.|<(  iG’(^)  —  G(r)'.  Therefore 

(7.5)  lim  —  ^11  =  0 

A’— oc 

From  AGy  =  (L(Gy),Gy)  and 

=  AP,^gy  =  (L{^B^G(t'^)y),  f  )  -  G(l'',))y 

^  k=l  fc=l 


it  follows  immediately  that  (because  Ylk=o  ^k^i^k)  —  ^i^G  and  P/^G  — >  G  with  respect 
to  the  graph  norm) 

(7.6)  lim  WA^'g^'  -Ag\\  =  o. 

S  -—oc 
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The  representation  (4.10)  together  with  Theorem  4.1, a)  shows  that  for  to  =  4r 

(7.7)  M-''’5-'^(<o)||  <  for  all  A', 

where  the  constant  M  does  not  depend  on  A^ 

Let  T^{t),  t  >  0,  denote  the  semigroup  generated  by  .  By  the  variation  of 

constants  formula  we  have 

(7.8)  T^’{t)  =  S^{t)-  f  S^{t-s)g^C^T^{s)ds,  t>0. 

Jo 

By  GronwaH’s  inequality  and  (7.5)  we  obtain 

(7.9)  \\T^’m\  <  const.ll^ll,  0  <  t  <  to,  iV  =  1, 2 - 

This  together  with  (7.G)  and  (7.7)  gives 

M-"’r'^'(#o)-’ll  <  A/l|cl|  +  const.||A"'a-"’||  lU-ll  <  M\\  =  l 

for  all  N,  where  M  is  not  dependent  on  N.  From  [20;p.56],  we  see 

(7.10)  = -g^'C^')y^dX.  t>2to. 

where  f  =  f]  ufaUfa  and 

fi  =  {A  1  ImA  =  ReA  <  /?}. 

f2  =  {A  I  jImA!  <  ReA  =  /?}. 

f3  =  {A|Uf,}. 

with  a  =  (l/to)ln2A7,  d  =  i  +  sup^  ||  (")  the  constant  defined  in  Lemma  2.4). 

Moreover  (see  [20:Thm.  4.4])  we  have 

(7.11)  (^^'  -^^C-^))"'||  <  const.jlmAI 

for  A  G  So,/?  =  {A  I  |ImA|  >  ReA  <  /?j,  where  ‘“const.”  does  not  depend  on 

N. 

For  z  £  Z  we  consider  the  equation 

(7.12)  yi-{A^  -g^c^)y^  =  P^z. 

Let  Ao  €  pi  A)  and  put 


Tx  =  /  +  (Ao/-^)-*({A-  Ao)/  +  ClC). 
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It  is  easy  to  see  that  A  £  p{A  —  QC)  if  and  only  if  is  continuously  invertible.  If  this  is 
the  case  then 

(7.13)  (a/  -{A-  goy'  =  -  -ArK 

Since  Ao  £  )  for  N  sufficiently  large,  we  analogously  have 

(7.14)  (a/-(>^^’ -a^’c^’))~'  =(r^r\Xoi  -A^’)-'^ 

for  A  €  p{A^  —  G^C^')  and  N  sufficiently  large,  where 

jrf  =  /  +  (Ao7  -  A^y'  ((A  -  Ao)/  +  G^’C^ . 

Using  I  =  Tx  —  iX(jl  —  A)~^  ((A  —  Ao)/  +  GC),  A  6  p{A  —  QC).  we  see  that  equation  (7.14) 
is  equivalent  to 

(7.15)  Tyi  +  s^')z^’  =  (Ao/  -  A^'r^p-^'z 


where 


(7.1G) 


sy  =  py  (( A  -  Ao)((Ao/  -  A^'r^P^'  -  (Ao/  -  A)-^) 

+  (Ao/  -  A’^'r^G'^'c^'  -  (Ao/  -  Aygc). 


Since  (Aq/  —  A)  Ms  compact,  we  have 

|lP''  (Ao/ -(Ao/-^)-'ll  ^0  a.s.V-.cc 

and  therefore  by  Lemma  B.l  also  [[(Aq/  —  A'’^  )~^  P^  —  (Aq/  —  A)~^  ||  — >  0.  Using  this  and 
(7.5)  we  obtain  from  (7.1G)  that 


li  0  as  A’  -+  oo 

uniformly  for  A  in  compact  subsets  of  p{A  —  GC).  Therefore  for  any  compact  subset  K  of 
p{A  —  GC)  there  exists  an  A'o  such  that  +  Cj^))  exists  for  A  £  K  and  N  >  Nq, 

i.e.,  K  C  piA^’  —  G^C^  )  for  A"  >  Aq.  Equations  (7.14)  and  (7.15)  imply 

OC 

for  e  K,  N  >  Ao. 

/fc=0 


This  representation  shows  that 

Jim  ||(ri:'r’-;r;'||  =  o 

A  —  OC 
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uniformly  on  compact  subsets  of  p{A  —  QC)  which  immediately  implies 
(7.17)  lim  \\(\I  ~{A^’ (XI -{A-gC))~'z\\^0 

N  —  ac 

uniformly  on  compact  subsets  of  p{A  —  gC).  The  last  result  implies  (by  the  same  arguments 
as  those  given  in  the  proof  of  Theorem  5.1)  that  we  can  replace  /?  in  the  definition  of  F  by 
any  uj  >  uj2  —  sup{ReA  |  A  €  cr{A  —  QC)]  and  still  have  the  representation  (7.10)  for  N 
sufficiently  large.  As  in  the  proof  of  Theorem  5.1  (but  using  (7.11)  and  (7,17)  instead  of 
Theorem  4.1)  we  obtain 

l|r^(t)|l  <  const.e‘"',  t>2to, 
which  together  with  (7.9)  (for  t  G  [0, 3to])  finally  gives 

(7.1S)  l|F'’^’(t)ll  <  const.e*^',  t  >  0, 

for  N  sufficiently  large.  | 

Corresponding  to  problem  (7.1)  we  consider  the  following  approximating  problems: 


Given  z  ^  Z  minimize 

(7.19)  J-"(t/;r)=  /  (lC-^’c^'(0|' +  lit(t)n  dt 

Jo 

subject  to 

(7.20)  z’^'(i)  =  S^'{t)P^'z  +  /  S^\t  -  s)Bu{s)ds,  i  >  0. 

Jo 

THEORE.M  7.4.  Suppose  that  (A,B)  is  stabilizablc  and  (C^A)  is  detectable.  Then  for  X 
sufficiently  large  the  algebraic  Riccati  equation  associated  with  (7.19). 

(7.21)  {A^'yu  +  UA^'  -  +  (C^'ye^’  =  0. 


has  a  unique  self-adjoint,  non-negative  solution  11'''.  There  exist  positive  constants  Mi, 
M2  and  u,'  such  that 

line’ll  <  Ml  and  <  Mje"'*'',  t  >  0, 

for  X  sufficiently  large. 

Proof:  The  result  on  existence  and  uniqueness  of  a  solution  11'^  of  (7.21)  for  X  sufficiently 
large  follows  immediately  from  Theorem  7.3  and  Theorem  7.2  applied  for  v4''^ .  B^  and  . 

Let  IC^  €  C(Z‘'^  iR”’)  and  G  £(R'’,2^  )  be  defined  as  in  the  proof  of  Theorem  7.3. 
Then  for  all  z  £  Z^ 


{U^'z,z)=  inf  J^'(u,z)<J^(-K^'z^\-),z) 
u6/,^ 


3G 


where  k  is  not  dependent  on  N.  Since  11^  is  non-negative,  this  implies 

<  K  for  all  N. 

From  equation  (7.21)  we  immediately  get 

-f  +  (C^’)T'^'  =  0. 

Let  2^(t)  =  exp((-4'''  —  B  )*n^’)f)2,  f  >  0.  Then  a  short  calculation  gives 

j^{z^\t).u^'z^'{t))  +  |(s^)*n^’2^(t)p  -H  |c^2^(t)l'  =  0,  t  >  0, 

or 


+  J^(j{B^'rU-’^'z^{s)\^  -h  lc^2^(s)f)  ds 

=  (n^’z,^)  <  ^llrlp.  f>0. 


As  a  consequence 


(7.22) 


-H  |C"’2'''(s)i')  ds  <  k||z|12,  t  >  0. 

By  the  variation  of  constants  formula  we  get 

=  r-^'(/)z  -f  [  T-^'u  -  s){g'^'C‘’^'  -  B^’{B-'^'yu-''')z^^\s)ds.  t  > 

Jo 

where  (t)  =  exp((v4-^  -  )f),  t  >0.  Taking  norms  and  observing  (7.5).  (7. IS)  and 

(7.22)  we  have 


0. 


ll--''(()ii  <  .i/<-'(ii--||  +  sup  lie'' II  +  /'ic'''.-'''(..)i* 

+  I|P|| /'|(e'Yn'.-'uoi*) 

<  .u,-'(it.-||  +  sup|ie''||('/»(jf“|r''.-''(.<)l"*)’'- 

+  liei|(’'"(/'|(g''rn''.-''(5)f  *)■'") 

<  .A/f“'(i  +  (‘/2(sup  lies'll  +  K''=iie||))ii.-|| 

for  all  f  >  0  and  all  .V.  Finally  we  get 

/“||.-''(()f<i(  <consl.||.-l|, 

Jo 

where  “const.”  is  not  dependent  on  N.  Now  the  result  follows  from  a  theorem  by  Datko 
(see  the  version  given  in  [22;Thm.  6.2]).  | 

The  following  theorem  is  an  immediate  consequence  of  Theorem  7.4  and  results  given  in 
[7;Thm.  6.9]  and  [4:Thm.  2.2].  These  results  eissume  that  the  adjoint  semigroups  converge 
strongly  as  is  the  case  for  our  scheme. 
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Theorem  7.5.  Suppose  that  (A.B)  is  stabilizable  and  (C^A)  is  detectable.  Then  the 
unique  solution  11'’''  of  (7.21)  converges  strongly  to  the  unique  solution  11  of  (7.2).  The 
semigroups  generated  by  A  —  are  exponentially  stable  uniformly  with  respect 

to  N  >  No  for  some  Nq  >  0.  Moreover,  for  any  e  >  0  there  exists  an  Nie)  such  that 

J[-B-U^'P‘^’z{-),  z)  <  J{u\  z)  +  ellrip 
for  allze  Z  and  N  >  N{e),  where  u°  =  -B*Ilz{). 
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Appendix  A 


In  this  appendix  we  state  and  prove  some  estimates  involving  the  Pade-approximant 

l  +  r/2  • 

Lemma  A.l.  For  m  >  1  and  |z|  <  2m 


e 


VH-z/2m/  1 


< 


(( 


l  +  |z|/27n 
l-|i|/2m 


_L_|,i3l±l£t/2m 
8m^  '  I  1  — |j|/2m 


for  Re  r  >  0. 


for  Rer  <  0, 


Proof:  Let  p=  z/m.  Then 


(A.l) 


For  \p\  <  2 


and  therefore 


(A.2) 


l-p/2 

l  +  />/2 


)' 


V  1  +  0/2/^  VI +  0/2/ 


m  — 1 


1  -p/2\'"-*-i 


l+P/2^,t^  Vl+^/2. 


l-p/2 
l  +  p/2 


jfc=t3 


=  -P_ 


1-P/2|^  1 


l  +  p/2  I  -  12 


k=A 


<gW(i+i;(2)  I'’' )  =  s''’' 


t-i 


For  Re  p  <  0  we  have 


m  — I 


l  -  p/2\"'-'=-i 


m  —  1 


it=0 


,-fcRep/l±H/r\"’ 

Vl-|p|/2>' 


-fc-1 


Jt=0 

^  -mRe^//l  +  IH/2\"*  \/1  +  1p1/2  \-l 

((rrwTs)  -")(r=w75-')  ■ 


This  together  with  (A.l)  and  (A.2)  gives 

_  /I  -  ^/2m\^|  1  /l^iy  1  +  kl/2m  R^,l  -  |^|/2m  //I  +  |-|/2m\^  _  \ 

vl  +  r/2m/  l~S\m/  1  -  b|/2m  Ul/m  \\l-|r|/2m/  / 


^1/2 
Nv.- 


1  +  lr|/2m\’" 
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for  l^l  <  2m  and  Re  z  <  0.  In  case  Re  p  >  0  we  get 

^=0  *  '  r  /  "  ik=0 

Therefore  using  again  (A.l)  and  (A. 2)  we  get 

T  —  z/2m\”'\  1 


-I  _  < 

\1  +  zl2m)  ~ 


l_e-Rep  - 
13  1  +  \z\/2m 


<  m. 


8m2'‘'  \-\z\l2m 


for  \z\  <  2m  and  Re^r  >  0.  | 


—  ^i.(.'rX/2A')  •  T'hen  for  jAj  <  p  and  N  sufficiently  large 


Lemma  A. 2.  Let  at  = 

(N  >  rp/2) 

for  Re  A  >0, 

'  ^pV(l  +  g)e'->((l±g^)''-l)  forReA<0, 

k  =  l,...,N. 

Proof:  Let  Ca’,a-  =  Xkr/N.  Then 


a  jt  <  { 


Ok  = 


Vi  +  CA-.ii72H  I’  . 


Qk  < 


For  |A(  <  p  we  have  |Ca’,*.  |  ^  pkr/N  and  |CA\ii  |  <  2A’  for  N  >  rp/2.  Lemma  A.l  for  r  =  Ca',*.- 
implies 

- 1)  to'  <  0- 

1  forReA>0, 

for  |A|  <  p  and  K  >  rp/2.  The  estimates  of  the  lemma  are  an  immediate  consequence.  | 
Corollary  A. 3.  For  any  compact  set  K  C  C  there  exists  a  constant  c  =  c{K]  such  that 

sup  |c>  (^)  —  j  <  c( A')(  — for  all  X . 

-r<0<Q  A 

Proof:  Since  ^k  is  the  first  order  interpolating  spline  for  .  we  have  (see 

[23;  Thm.  2.6]) 


-r<9<0 


k=0 


-r<9<0 


Then  the  result  follows  from  the  estimate 

l-^A/2-^^^^ 


Xil"  /  ^  ~  r  A/2A  ^  k 


<  max  c'''*  —  (r 
k=l . A'  '1 


+  rA/2A^)‘  <const.(^)’ 


for  all  A  €  A'  and  A’'  =  1, 2, . . .  ,  which  is  an  immediate  consequence  of  Lemma  A. 2.  | 
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Corollary  A. 4.  For  any  compact  set  K  C  C  there  exists  a  constant  c  =  c(A')  such  that 
lA-^’(A)  -  A(A)|  <  c{K){^f  for  A  e  K  and  A  =  1,2, . . .  . 

PROOF:  By  definition  of  A^’(A)  and  A(A)  we  have 

|A''(A)-A(A)|<7  sup  ie«(«)-e-'«|, 


where  7  is  the  constant  defined  in  Lemma  2.4.  Then  the  result  follows  immediately  from 
Corollarj'  A. 3.  I 


Lemma  A. 5.  There  exist  positive  constants  A'o  and  ao  such  that  for  0  <  j  <  N,  N  = 
1,2, . . . ,  and  A  €  Da.u; 


1  -  rX/2N  ^ 
l  +  rX/2N 


<  KqC 


‘|ImA| 


provided  a  >  oq.  In  fact  A'o  =  and  oq  =  ^(In  -  l)  is  a  possible  choice  for  A'o 

and  Oq. 


Proof:  Let  Oo  =  ^(1  -  XVe  consider  the  following  three  cases  separately. 

Case  1:  0  <  ReA  <  w.  A  £ 

In  this  case  the  result  follows  immediately  from 


1  -  rA/2.V 


1  +  rA/2A' 


<  1  <  €' 


^llmAI  <c‘-‘'e-“^|ImA|. 


Case  2:  — qq  <  Re  A  <  0,  A  €  Sa.u-- 

We  put  X  =  Q  -hid-  The  inequality  a  >  —ao  implies  1  +  (r/2A)Q  >  0  and  therefore 


1  -  rA/2.V 


1  +  rA/2.V 


(1  -  ra/2Ny  +  (r^/2A’^V^^  <  1  " 


(1  +rQ/2A')2  +  (r5/2A')2 


1  +  ro  /2.V 


For  the  function  /(a)  =  ( l  +  ro/2A  )  -2N/r  <  o  <  0,  we  have /(O)  =  1  and  /'(a)  >  0 

for  — Oo  <  a  <  0.  Therefore  /(a)  <  1  on  —00  <  Q  <  0  and 


1  -  rA/2A’ 

2 

1  -  rA/2A 

1  +  rA/2A' 

H-rA/2A' 

N 


<  e 


-6Re  A 


<e-“^|ImA|,  j  =  l.....A. 


Case  3:  ReA  <  — Qo,  A  G  Ha..*- 
Let  A  =  —a  +  t/3,  o  >00-  Then 


(A.3) 


l-rX/2A' 


1  +  rA/2.V 


iN 


1  - 


1  +  rA/2.V  I 


<  1  + 


ll  +  rA/2.V| 


(  2  \^  (  1  4.Y2\ 


VI 
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We  claim  that 


(A.4) 


— <  ba  for  A  e  a  >  qq,  A’’  =  1, 2, . . . , 

rl^l 


provided  a  is  sufficiently  large.  Then  we  get  from  (A. 3) 


1  -  rX/2N  ^  ^  1  -  rX/2N 
l  +  rX/2N  -  l+rX/2N 


<(l  +  ^)'"<e‘“<e-“‘|ImA| 


for  j  =  1, . . . ,  iV  and  a  sufficiently  large. 

We  now  prove  (A.4).  Since  \(3\  >  the  estimate  (A.4)  is  certainly  satisfied  if 

/r  <  6Qe*“e“*’,  a  >  Oq.  Since  the  function  are*  is  increasing  on  a:  >  0,  (A.4)  is  satisfied 
if 


(A.5) 


<  6ao£*“‘’e“\ 


The  definition  of  oq  shows  that  (A. 5)  is  equiralent  to 

Using  ae~*  <  e”^  for  x  >  0.  we  see  that  the  last  inequality  is  satisfied  for  all  .V  provided 
e"**  >  r(6(6  —  r)e)  ^ ,  i.e., 
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Appendix  B 

In  this  appendix  we  prove  various  rate  estimates  for  the  convergence  of  the  resolvents 
corresponding  to  the  generators  of  the  approximating  semigroups  We  choose  k  G  R 

such  that  0  G  p{A  —  kT).  Then  Theorem  4.1,b)  shows  that  also  0  G  p(A^  —  kT‘'^  )  for  N 
sufficiently  large,  N  >  Ni .  We  put 

A  =  A-kT  and  A^  =  A^  -  kT^  ,  N  >  Ni. 


Furthermore  we  introduce  the  notation 


f  /-O 

L{if)  =  {Ao-Ki„)^io)  +  Y.^M^j)+  /  A{e)^{e)dd, 

j=i 


A'^'(A)  =  A/„  -  A(A)  =  A/„  - 


.N 


A 


The  following  formula  will  be  of  fundamental  importance: 

P^  (A7  -  A)-^  -  (A/  -  A^r^P^ 

=  (XI  -  A^')-'^  (^P^{A  -  kT)  -  (A^  -  kT^)P^')(A/  -  A}-^ 

=  (XI  -  A"'’ )-^ A^ P^' (a-^ P^’ A  -  I^XI  -  A)-'^ 

=  ((A/  -  >l-“^')-i(A7  -  kT^’)  -  I^P^{A-^P^'A  -  7)(A7  -  A)'^ 

for  A  G  p{A)  and  N  >  N\.  Here  we  also  have  used  A^' z  =  AL^  z  for  c  G  Z  '.  The 
constants  a',  b  and  a  are  chosen  as  in  Section  6. 

Lemma  B.l.  Let  a  =  1  if  (H)  is  satisfied,  otherwise  a  =  1/2. 
aj  There  exist  positive  constants  c  and  A’l  such  that 

WP-'  iM  -  A}-'  -  (A/  -  <  c(^)“|Im  A|< 

for  aJJ  X  G  Tg  and  ,V  >  A'] . 

bj  For  any  compact  set  K  C  p{  A)  there  exist  positive  constants  c  and  A'l  such  that 

11P''(A/  -  Al)->  -  (XI  -  .4^’)->pA'||  <  c(^)“ 
for  all  X  G  K  and  N  >  Aj . 

PROOF:  For  (p  S  (i.e.,  (:,^(0),v5)  €  dom>l)  we  get 

pX'A(^(0M  =  P'-'(L(v),^)  =  (L('^).x^') 


with  (see  Lemma  2.1, a)) 


A' 


k=l 
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From  (1.5)  we  obtain 

(B.2)  A-^P^’Ai^(0\^)  = 

where  (observe  A'’^  (0)  =  A(0)) 

'A''W  =  'A''’(0)+ /Vm*, 

<^"(0)  =  A(0)-'(£(4»)  +  i( /\'*(r)rfr)). 

It  is  obvious  that  x^{0)  :=  Jg\^’{T)dT  =  Y,k=i  BkiB)x’^{t^)  with 

X^’iti^)  =  <p(0)  -  k=l,...,N. 

This  shows  x'^  =  (^(0)  -  ,  where  (p^  =  Bj^’tp{t^).  Therefore 

T/.^(0)  =  A(0)-'i(<^  -p^  +  <pi0))  =  A(0)-'L(<^  -  p^)  -  p{0), 

=  A(0)~'L(c^  ^(0)  +  9(0)  -  9^^'  =  A(0)-'i(s?  -  9^^')  -  9""'. 

This  together  with  (B.2)  gives 

(B,3)  (.4-‘P«i  -  =  -(A(O)-’i(v  -  A(0)-'i(^;  -  9''')  +  9  -  9-'''). 

If  (H)  is  satisfied  we  use 

(B.4)  |i(9-9''’)l  =  I /°.'l(«)(9(«)-9'''W)rf»|  <l|.4|U.||9-9'"''lli.>, 

otherwise  (7  =  |,4„  -  k/„|  +  |.4;|  +  |.4(»)|  tK) 

(B.5)  Ii(9-9^’)l<7ll9-9''lk«. 

In  case  of  (B.4)  we  use  ([23;  Thm.  2.4]) 

||9  -9^IIl»  <  const. 
in  case  of  (B.5)  ([23;  Exercise  (2.10)]) 

||9  -  9^’IU«  <  const. (^)^ ''^119 II , 

in  order  to  get  from  (B.3) 

||(.4-'P-'U  -  /)(9(0),9)ll  <  const.(^)“||9||„. 
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or,  equivalently, 


(B.6)  -  711c(d„„. AZ)  <  )° 

for  N  =  1,2,...  . 

From  Theorem  4.1, a)  and  the  estimate  |A|  <  const. |ImA|  for  A  €  ^a,u/  obtain 
(B.7)  ((A/-^^)-^(A/-^£7'^)-/)p^|  <const.|ImA|2 

for  all  A  G  Sa.u;  and  N  >  Aq.  Using  A{XI  —  A)~^  =  A(A/  —  -4)”*  —  /  we  see  that  also 

((^(A/  —  <  const.(Im  Ap||r|| 

for  all  X  £  Ha, ui  and  z  £  Z.  Therefore 

(B.S)  ||(A/ -  ^)"^l|c(Z,dom^)  <  const.jlmAp  for  A  6 

Now  part  a)  of  the  lemma  follows  from  (B.l)  and  (B.6)-(B.S). 

In  order  to  prove  part  b)  we  use  Theorem  4.1,b)  in  order  to  see  that  the  right-hand  side 
in  the  estimate  (B.7)  can  be  replaced  by  a  constant  as  long  as  A  £  K  and  N  >  AV  For 
the  estimate  (B.S)  this  is  obvious.  | 

Lemm.a  B.2.  Let  /?  =  2  if  (H)  is  satisfied,  otherwise  0  =  3/2. 

a)  There  exist  positive  constants  c  and  Nj  such  that 

|1P-^'(A/  -  Aj-K-  -  (A/  -  A^'r^P^'z]]  <  c(^)^|Im  Apl|r|ldon.  .4 

for  all  A  £  Ha.u.-‘  ^  €  dom^l  and  N  >  Ni.  If  z  =  {rj,^)  with  — r,0;  R"  )  then  the 

right-hand  side  of  the  estimate  is  c(-^)  jImAp((i7|  -I-  (((,?|(vvi.<=c). 

b)  For  any  compact  set  K  C  p{A)  there  exist  positive  constants  c  and  A’l  such  that 

for  all  A  £  K,  z  £  dom^l  and  N  >  A’j.  If  z  =  (»7,<^)  with  <p  £  U''^’°°(— r,0;  R”)  then  the 
right-hand  side  of  the  estimate  is  c(-^)  (jr/l  +  |li^||w>  »)- 

Proof:  For  z  =  (ri,ip)  we  put  (t/’(0),V^)  =  (A/  —  A)~^{T],ip).  In  analogy  to  (B.3)  we  get 
(B.9)  {A-^P’^A  -  /)(i/’(0),  V’)  =  -(A(O)-'Z(V'  -  0^),  A(0)-'Z(^  _  V>^) 
where  r/’^'  =  Ylk=o  ^^sing  (4.1)  we  get  (compare  the  derivation  of  (4.2)) 

|A(A)-'1<7]\7  forA€3:o...,  A^  =  1,2 . 

t)|Al 
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Then 


!V'(0)!  £  +7  sup  ]/ 

0|A|  \  -r<$<0  J»  ^ 

<  const. ll^ll  for  A  €  Ea,w 
This  together  with  (4.8)  gives 

(B.IO)  |0(^)I  ^  const. Ijlm  A|  ||2j|  for  A  €  and  —r<6<0. 

If  (p  is  absolutely  continuous  we  have  by  definition  of  xp 
(B.ll)  ij)  =  Xxl'  —  pi  ^  =  Xxp  —  p. 

Observing  jA]  <  const. jlmAj  on  Eg,^-  we  get  from  (B.ll) 

(B.12)  1Iv'’I1l~  <  const-jlm  Ap|li:||  for  A  e  Eq,^. 

and 

(B.13)  IKiU‘=‘  <  const. lImAp(|r?l  +  j|^l|ivi,c»)  for  A  6  Ea..^. 

if  p  6  H’^’°°(— r,0;  R”),  resp. 

(B.14)  |[t/>||is  <  const. |ImA(^(|77|  +  llc^llivi.j)  for  A  6  Da.*-, 

if  r,0:  R").  (B.14)  specifically  implies 

(B.15)  <  const. lImA|^||r||doni>t  for  A  G  E^.;^., 

if  ?  €  dom>t.  Depending  whether  t'’  is  in  or  in  we  have  the  estimates 

(B.16)  ||j/>  -  <  const. (^)^||v'’||w*, 

(B.17)  ||v-t/’^IU«  <  const. (^)^||v>||vv2.-, 

(B.18)  <  const.(^)’''^||^||HS 

(see  [23;  Thm  2.5,  Thm  2.6  and  Exercise  (2.10)]).  If  (H)  is  satisfied  we  use  (B.16)  together 
with  (B.14),  otherwise  we  have  to  take  (B.17)  or  (B.18)  together  with  (B.13)  or  (B.14), 
respectively,  in  order  to  get  the  estimates  of  part  a)  from  (B.9)  and  (B.7). 

The  proof  for  part  b)  is  analogous  but  simpler.  As  in  the  proof  of  Lemma  B.l  one  has  to 
use  the  fact  that  the  right-hand  side  of  (B.7)  can  be  replaced  by  a  constant  for  A  £  K .  | 
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